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Summary 

This  project  was  centered  around  research  and  development  of  next  generation 
metamaterials  that  are  needed  in  radar  antenna  technologies.  By  directing  electromagnetic 
beams,  such  as  used  in  radar,  such  materials  can  enormously  improve  detection  and 
communication  in  aircraft,  ranging  from  communication  with  satellites  to  detection  of  enemy 
bodies. 

During  this  performance  period,  groundbreaking  results  were  obtained  in  metamaterials 
that  display  negative  refraction  properties.  New  metmaterial  designs  were  developed  and  new 
phenomena  arising  from  the  negative  refraction  properties  were  demonstrated.  Our  work  in  this 
area  has  led  to  some  major  results,  reported  in  10  completed  papers,  and  has  received  extensive 
attention  in  the  popular  press,  including  writeups  in  Science,  Physics  Today  and  APS  News. 

Negative  Refraction  and  Left-Handed  Electromagnetism 

Naturally  available  materials  that  are  transparent  to  electromagnetic  (EM)  waves  can  be 
characterized  by  a  positive  refractive  index  n  >  0  .  However,  as  noted  by  Veselago  in  1968, 
Maxwell’s  equations  allow  the  possibility  of  negative  permttivity  s  <  0  and  permeability  p  <  0 , 

and  correspondingly  negative  refractive  index  n  =  -J\  e  ||  p  \  <  0  .  Recently  Smith,  et.  al. 

demonstrated  experimentally  that  in  certain  metamaterials  microwaves  bend  negatively 
corresponding  to  n  <  0  . 

Our  work  can  be  classified  into  three  areas: 

Theoretical  Issues  in  Negative  Refraction  and  Imaging 

The  most  striking  property  of  EM  waves  in  negative  index  media  (NIM)  is  that  of  Left- 
Handed  Electromagnetism  (LHE),  since  the  electromagnetic  fields  E  and  H  ,  and  the  wave 
vector  k  ,  form  a  left-handed  triplet  set.  Consequently  the.  energy  flow  represented  by  the 
Poynting  vector  S  =  ExH  is  anti-parallel  to  the  wave  vector  k ,  so  that  S -k  <  0 .  For  conventional 
n  >  0  materials  E,H,k  form  a  right-handed  set  with  S-k>  0 . 

Refraction  of  energy  in  wave  groups  Since  all  physical  radiation  sources  emit  their 
radiation  in  the  form  of  wave  packets,  we  analyzed  refraction  of  electromagnetic  wave  packets  on 
passing  from  an  isotropic  positive  (PIM)  to  an  isotropic  NIM.  We  have  treated  several  exhaustive 
examples:  localized  wave  packets,  beams,  and  also  a  finite  number  of  plane  waves.  In  all  of  the 
cases  we  show  that  the  energy  and  momentum  refract  negatively.  We  show  that  in  all  cases  where 
the  wave  group  does  not  extend  to  infinity  in  the  direction  perpendicular  to  the  wave  vector,  the 
interference  pattern  also  refracts  negatively.  This  paper  was  published  Physical  Review  E,  69, 
026604  (2003).  (Several  simulations  and  movies  are  available  at  our  website 

http://sagar.phvsics.neu.edu). 

Near  and  far  field  imaging  by  n<0  and  e  <  0  materials  We  studied  analytically  and 
numerically  imaging  by  Left-Handed  Metamaterials.  For  the  far  field,  the  image  is  “perfect”  only 
for  n  =  -1 .  Away  from  this  value,  aberration  and  caustics  will  be  present  and  were  analyzed.  For 
near  field  imaging  with  Re(s)  = -l,Re(/0  = +1  >  the  consequences  of  departure  from  dielectric 
matching  and  losses  were  analyzed.  This  paper  has  been  published  in  Microwave  and  Optics 
Technology  Letters,  39,  282  (2003). 
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Figure  1  Negative  refraction  in  a  Left-handed  Metamaterial.  From  Left:  Metamaterial  constructed  at  NU. 
Schematic  of  refraction  setup.  (Middle)  Angular  and  Spectral  dependence  of  refracted  signal  showing  regions  of 
positive  and  negative  refraction.  (Right)  Negative  refraction  in  LHM  compared  with  positive  refraction  in  Teflon. 


Negative  Refraction  in  Left-Handed  Metamaterials 

We  fabricated  negative  index  metamaterials  by  interleaving  arrays  of  split  ring  resonators 
and  wire  strips.  Microwave  measurements  were  carried  out  in  1-D  X-band  waveguide  and  2-D 
waveguide  setups.  From  the  ID  waveguide  results  we  determined  the  crucial  material  parameters 
such  as  complex  refractive  index,  permeability  and  permittivity  ( n  =  ri  +  in",  fu  =  ju'  +  iju"  and 

£=s'  +  is")  required  to  design  materials  and  optimize  characteristics  for  microwave 
applications.  The  n'  determined  from  the  waveguide  transmission  parameters  is  found  to  vary 
from  -1.9  to  -0.3  in  the  passband  region  9.6-10.5GHz.  The  results  show  that  transmission  is 
optimized  for  n'  =  -1  and  low  «" .  One  of  our  significant  contributions  in  the  field  of  left  handed 
electromagnetism  has  been  in  proving  its  physical  reality  in  both  NIM  and  photonic  crystals,  a 
topic  which  has  been  recently  controversial. 

Microwave  Photonic  Crystals 

Our  search  for  novel  materials  that  display  negative  refractive  index  have  led  to  the 
observation  of  negative  refraction  in  a  new  class  of  materials  -  metallic  photonic  crystals.  We 
have  demonstrated  negative  refraction  of  microwaves  in  metallic  photonic  crystals  fabricated 
using  cylindrical  copper  rods  arranged  on  a  periodic  lattice.  With  parallel  theory  and  numerical 
simulations  we  showed  that  the  negative  refraction  observed  corresponds  to  left-handed 
electromagnetism  and  arises  due  to  the  dispersion  characteristics  of  waves  in  the  periodic 
medium.  The  ease  of  fabrication,  control  of  design  parameters,  scalability  and  very  low 
attenuation  in  the  metallic  photonic  crystal  (100  times  smaller  than  that  in  composite 


Figure  2  (Top  left)  Representative  wave  field  scans  of  the  refracted  microwave  transmission  intensity  by  photonic 
crystal  prism,  demonstrating  (left)  positive  refraction  at  6.62  GHz,  and  negative  refraction  at  9.77  GHz.  The 
experimental  results  are  in  good  agreement  with  band  structure  calculations  and  numerical  simulations  of  wave 
refraction.  (Right)  Simulations  of  wave  refraction  showing  (a)  positive  refraction  at  7  GHz  and  (b)  negative 
refraction  at  9  GHz. 
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metamaterial)  pave  the  way  for  several  interesting  applications  from  microwave  to  optical 
frequencies.  These  results  was  published  in  Physical  Review  Letters,  92,  127401  (2004). 

Tailor-made  negative  refractive  indices  of  photonic  crystals 

Negative  refraction  and  left-handed  electromagnetism  in  a  metallic  photonic  crystal  (PC) 
wedge  were  demonstrated  in  free  space  for  both  transverse  magnetic  and  electric  mode 
propagation.  The  experimental  results  are  in  excellent  agreement  with  numerical  calculations 
based  on  the  band  structure  with  no  fit  parameters  used  in  modeling.  The  results  demonstrate 
precision  control  on  the  design  and  fabrication  of  the  PC  to  achieve  tailor-made  refractive  indices 
between  -0.6  to  +1.  These  results  were  published  in  Applied  Physics  Letters,  85,  1858  (2004). 

Imaging  by  a  flat  lens  due  to  negative  refraction 

The  positive  index  of  refraction  in  conventional  optical  lenses  requires  them  to  have 
curved  surfaces  to  form  an  image.  However,  negative  index  of  refraction  allows  a  flat  slab  of  a 
material  to  behave  as  a  lens  and  focus  electromagnetic  waves  to  produce  a  real  3-D  image  [see 
our  website  http://sagar.phvsics.neu.edu  for  an  introduction  to  this  idea].  We  have  demonstrated 
this  unique  feature  of  imaging  by  a  flat  lens,  using  the  phenomenon  of  negative  refraction  in  a 
photonic  crystalline  material. 

The  key  advance  that  enabled  this  observation  is 
the  design  of  a  suitable  photonic  crystal  to  achieve 
negative  refraction  over  a  wide  range  of  angles.  While 
negative  refraction  was  demonstrated  at  microwave 
frequencies  in  a  quasi-homogeneous  metamaterial, 
imaging  using  a  flat  lens  of  such  material  has  not  been 
successful  due  to  the  large  dissipation  and  anisotropy  in 
these  materials  and  the  lack  of  proper  design  tools.  This 
work  capitalizes  on  our  earlier  demonstration  of  negative 
refraction  in  photonic  crystals. 

Figure  3  Image  of  Microwave  Point  Source 

The  accompanying  figure  shows  the  image  of  a  microwave  point  source  of  frequency  9.3 
GHz  (wavelength  3.22  cm)  placed  at  a  distance  of  2.25  cm  from  a  2-D  flat  lens  made  of  a 
photonic  crystal  fabricated  from  an  array  of  cylindrical  alumina  rods.  On  the  far  side,  an  image  of 
excellent  quality  is  seen  at  a  distance  2.75  cm.  This  image  is  observed  only  in  a  narrow  frequency 
range  between  9.0  -  9.4  GHz,  with  the  best  focus  achieved  at  9.3  GHz.  Outside  this  narrow  band, 
at  all  other  frequencies  between  2-12  GHz,  a  single  focus  point  is  not  observed.  These 
observations  can  be  completely  understood  by  examining  the  band  structure  of  the  photonic 
crystal  and  corresponding  equi-frequency  surfaces  (EFS). 

Conventional  optical  lenses  have  a  single  optical  axis,  limited  aperture  and  cannot  focus 
light  onto  an  area  smaller  than  a  square  wavelength  [1],  In  contrast  the  present  flat  lens  does  not 
have  a  unique  optical  axis  and  is  not  restricted  by  the  aperture  size.  We  have  demonstrated  both 
these  features  by  moving  the  object  by  4  cm:  the  image  moves  by  a  corresponding  amount  in  the 
same  direction  (see  accompanying  figure).  Note  for  the  subwavelength  source  we  have  observed 
an  image  of  similar  size.  The  image  formed  by  the  flat  lens  due  to  negative  refraction  is  a  real  3D 
image. 
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Figure  4  A  flat  lens  without  a  uniqe  optical  axis 
Shows  an  image  moves  the  object  by  4cm 

The  unique  properties  of  the  flat  lens  reported  here  lead  to  entirely  new  perspectives  on 
imaging  and  potentially  new  applications.  A  particular  advantage  of  the  photonic  crystalline 
material  is  its  scalability  to  sub-micron  dimensions  ensuring  several  possible  applications  from 
microwave  to  optical  frequencies.  The  paper  describing  this  work  was  published  in  Nature,  426, 
404  (2003).  This  work  was  subsequently  cited  by  Science  as  one  of  the  breakthroughs  of  the  year 
2003. 


Focusing  by  Planoconcave  Lens  Using  Negative  Refraction 

We  demonstrated  focusing  of  a  plane  microwave  by  a  planoconcave  lens  fabricated  from 
a  photonic  crystal  having  a  negative  refractive  index  and  left-handed  electromagnetic  properties. 
An  inverse  experiment,  in  which  a  plane  wave  is  produced  from  a  source  placed  at  the  focal  point 
of  the  lens,  is  also  reported.  A  frequency-dependent  negative  refractive  index,  n(w)<  0  is  obtained 
for  the  lens  from  the  experimental  data  which  match  well  with  that  determined  from  band 
structure  calculations.  This  work  was  published  in  Applied  Physics  Letters. 

Slow  Microwaves  in  Left-handed  Materials 

Remarkably  slow  propagation  of  microwaves  in  two  different  classes  of  left-handed 
materials  (LHM’s)  was  demonstrated  from  microwave-pulse  and  continuous-wave  transmission 
measurements.  Microwave  dispersion  in  a  composite  LHM  made  of  split-ring  resonators  and  wire 
strips  reveals  group  velocity  vg~c/50,  where  c  is  the  free-space  light  velocity.  Photonic  crystals 
(PhC’s)  made  of  dielectric  A1203  rods  reveal  vg~c/10.  Group  delay  dispersion  of  both  the 
composite  LHM  and  PhC’s  determined  from  the  experiment  is  in  complete  agreement  with  that 
obtained  from  theory.  The  slow  group  velocities  are  quantitatively  described  by  the  strong 
dispersion  observed  in  these  materials.  This  work  was  published  in  Physical  Review  B. 


Flat  lens  without  optical  axis:  Theory  of  imaging 

We  derived  a  general  theory  for  imaging  by  a  flat  lens  without  optical  axis.  We  showed 
that  the  condition  for  imaging  requires  a  material  having  elliptic  dispersion  relations  with 
negative  group  refraction.  This  medium  is  characterized  by  two  intrinsic  parameters  a  and  k. 
Imaging  can  be  achieved  with  both  negative  and  positive  wave  vector  refraction  if  o  is  a  positive 
constant.  The  Veselago-Pendry  lens  is  a  special  case  with  a  =  1  and  K  =  0.  A  general  law  of 
refraction  for  anisotropic  media  is  revealed.  Realizations  of  the  imaging  conditions  using 
anisotropic  media  and  inhomogeneous  media,  particularly  photonic  crystals,  are  discussed. 
Numerical  examples  of  imaging  and  requirements  for  sub-wavelength  imaging  are  also  presented. 
This  work  was  published  in  Optics  Express. 
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Flat  lens  without  optical  axis:  Theory  of 

imaging 
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Abstract:  We  derive  a  general  theory  for  imaging  by  a  flat  lens  without 
optical  axis.  We  show  that  the  condition  for  imaging  requires  a  material 
having  elliptic  dispersion  relations  with  negative  group  refraction .  This 
medium  is  characterized  by  two  intrinsic  parameters  a  and  k.  Imaging 
can  be  achieved  with  both  negative  and  positive  wave  vector  refraction  if 
cr  is  a  positive  constant.  The  Veselago-Pendry  lens  is  a  special  case  with 
a  =  1  and  K  =  0.  A  general  law  of  refraction  for  anisotropic  media  is 
revealed.  Realizations  of  the  imaging  conditions  using  anisotropic  media 
and  inhomogeneous  media,  particularly  photonic  crystals,  are  discussed. 
Numerical  examples  of  imaging  and  requirements  for  sub-wavelength 
imaging  are  also  presented. 

©  2005  Optical  Society  of  America 

OCIS  codes:  (110.2990)  Image  formation  theory;  (260.2030)  Dispersion;  (100.6640)  Super¬ 
resolution;  (080.3620)  Lens  design 
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Since  antiquity,  the  positive  index  of  refraction  of  conventional  materials  has  required  the  use 
of  curved  surfaces  to  focus  light.  However  there  has  been  a  continuing  quest  for  lenses  with 
flat  surfaces  as  they  confer  a  variety  of  advantages.  Notable  examples  are  the  Fresnel  lens  and 
Maxwell’s  fish  eye  lens  [1, 2].  The  Fresnel  flat  lens  uses  a  gradient-index  n(x)  material  (GRIN) 
[3,  4]  and  hence  must  possess  an  optical  axis,  i.e.  is  not  translationally  invariant  along  the 
surface  (x-axis). 

The  concept  of  negative  refraction  has  led  to  new  fundamental  approaches  as  well  as  appli¬ 
cations  in  optics  [5].  In  1968  Veselago  [6]  pointed  out  that  for  a  material  with  refractive  index 
n  =  —  1  a  flat  surface  would  focus  light.  A  decade  later,  Silin  [7]  discussed  a  more  general  case 
and  obtained  the  lens  equation  of  a  flat  slab  with  negative  elliptic  dispersion.  Pendry’s  recent 
analysis  [8]  demonstrating  the  possibility  of  sub-wavelength  resolution  with  a  flat  slab  of  such 
materials,  as  well  as  the  experimental  realization  of  the  so-called  left-handed  materials  using 
composite  media  [9]  and  photonic  crystals  (PhCs)  [10,  11]  led  to  renewed  interest  in  the  unique 
electromagnetic  properties  of  these  artificial  materials.  The  unique  property  of  these  flat  lenses 
is  the  lack  of  optical  axis.  This  type  of  flat  lens  can  be  realized  in  a  PhC  using  negative  refraction 
[12]  and  has  been  demonstrated  in  microwave  experiment  [13].  There  is  however  no  complete 
theory  of  imaging  by  a  flat  lens  which  properly  describes  the  various  features  observed  in  the 
experiments  and  in  numerical  simulations. 

In  this  paper,  we  present  a  general  theory  of  imaging  by  a  flat  lens  without  optical  axis, 
resulting  in  a  proposal  for  a  new  material  which  is  translationally  invariant  along  the  surface 
and  has  anisotropic  refractive  index  n(6).  Defining  an  optical  phase  condition  for  imaging,  we 
show  that  the  condition  requires  a  material  having  Elliptic  dispersion  relations  with  Negative 
Group  Refraction  (ENGRM)  at  the  operating  frequency.  The  ENGRM  is  defined  by  two  pa¬ 
rameters:  the  anisotropy  parameter  <7  which  is  a  measure  of  the  ellipse  eccentricity,  and  a  phase 
factor  k  which  determines  the  center  of  the  ellipse.  Refraction  laws  for  the  wave  vector  and 
the  group  velocity  for  anisotropic  media  are  derived.  The  theory  shows  that  “perfect”  images 
can  be  obtained,  and  imaging  is  possible  with  negative  as  well  positive  refractive  indices .  The 
required  conditions  and  consequences  of  flat  perfect  and  imperfect  lenses  made  of  ENGRM  are 
derived  from  a  generalized  Fermat’s  principle.  Realization  of  flat  lens  using  PhCs  is  discussed. 
In  real  PhCs,  <7  is  itself  angle-dependent  in  most  cases,  leading  to  some  limitations  for  image 
formation.  Many  of  the  proposals  for  PhC  lenses  are  contained  in  the  present  theory.  The  the¬ 
ory  is  also  applied  to  recent  experiments  using  PhC  and  is  shown  to  successfully  describe  some 
intriguing  aspects  of  the  data.  Numerical  simulations  are  presented  that  describe  visual  details 
of  the  image  formation  in  PhC. 

We  first  consider  the  case  of  a  single  interface  which  is  along  the  x-axis  and  at  z  =  0.  The 
wave  vector  k  of  an  incident  plane  wave  from  vacuum  towards  the  interface  makes  an  angle  9 
with  the  z-axis.  Here  we  consider  the  case  that  for  any  incident  k,  there  is  only  one  refractive 
wave  vector  kr.  Continuity  along  the  interface  requires  that  krx  —  kx.  A  point  source  is  at  x  =  0 
and  z  =  —u  with  u  >  0.  We  now  consider  another  point  in  ENGRM  at  x  =  Ax  and  z  =  v  with 
v  >  0.  The  phase  difference  between  these  two  points  for  each  incident  k  is  O  =  Oz  +  O^with 
®z  =  kzu+  knv  and  <t>x  =  kxisx.  According  to  Fermat’s  principle,  the  formation  of  an  image 
would  require  that  the  total  phase  O  be  stationary.  This  restriction  can  be  relaxed  for  a  flat 
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interface.  Since  the  surface  is  flat,  any  point  of  a  finite  size  source  will  be  imaged  to  a  point  with 
the  same  ^-coordinate.  Thus  Ax  =  0  and  O*  =  0.  One  only  needs  to  consider  ®z.  A  generalized 
Fermat’s  principle  states  that  an  image  will  be  formed  if  the  phase  <f>z  is  stationary,  dOz/ dkz  =  0. 
The  lens  equation  for  single-interface  is 

u  =  av  (1) 

with  a  material  constant 

a  =  —dkrz/dkz.  (2) 

This  constant  tr  should  be  positive  and  independent  of  the  incident  angle  for  a  focus  with¬ 
out  aberration.  Thus  the  following  rule  for  the  wave  vector  refraction  at  the  interface  must  be 
obeyed 

krZ  ~  K  okz.  (3) 

Here  k  is  the  integration  constant  which  is  also  an  intrinsic  property  of  the  ENGRM.  Since  in 
the  vacuum  k2  +  k2  =  k$  —  a%/c2,  the  equi-frequency  surface  (EFS)  of  the  medium  is 

a~2(K-krz)2+k2x  =  (4/c2.  (4) 

Note  that  this  elliptic  dispersion  exists  only  at  the  operating  frequency  ft)o  [14].  In  the  neigh¬ 
borhood  of  this  frequency  <w  ~  too,  ( K  —  krz)2/n2  +  k2x/n 2  =  <o2/c2  with  cr  =  nz/nx  ~  a.  The 
deduction  of  Eq.  (3)  from  an  elliptic  EFS  requires  negative  group  refraction 

(kr  -  Ki)  ■  Vkrm  <  0  (5) 

together  with  the  causality  condition  d(o/dkrz  >  0  that  the  energy  must  flow  away  from  the 
interface  which  is  illustrated  in  Fig.  1.  The  ray  vector  [15]  which  represents  the  direction  of 
group  velocity  is  s  =  (— a  sin  ax  +  cos  az)/nxa  while  kr  —  fez  =%(sin  ax  —  cfcosaz)  with 
sin  a  =  nf1  sin0. 

The  Snell’s  law  for  wave  vector  loses  its  meaning  in  this  anisotropic  medium.  However  there 
is  a  law  for  group  refraction,  which  is  tan /3  =  —  crtana  =  — cr(n2  —  sin2 0)-1/2sin0.  At  the 
operating  frequency  cod  of  the  flat  lens,  nx=  1,  ct(coo)  =  a,  the  group  refraction  law  is  simply 

tan/3  =  —  crtan0.  (6) 

Thus  —a  can  be  regarded  as  the  effective  refractive  index.  Note  that  this  refraction  law  is 
very  general  and  is  valid  for  any  a  as  defined  by  Eq.  (2).  The  proof  is  the  following,  tan  j3  = 

=  - Here  we  used*«  =  , 

The  space  inside  the  medium  can  be  considered  to  be  “optically  stretched”  in  the  z-direction 
by  a  factor  of  a-1.  The  waves  inside  the  medium  are  represented  by  i/r(r')  =  exp(ikxx  —  ikzf  + 
iidz!)  with  f  =  Oz  and  K1  =  k/o.  Thus  in  the  case  that  k  0,  the  waves  at  this  frequency 
will  experience  a  periodicity  in  the  z-direction  with  length  scale  2k/\k\  since  the  waves  are 
modulated  by  exp(ticz).  This  periodicity  is  the  effective  periodicity  and  may  not  be  the  actual 
physical  periodicity  such  as  the  lattice  spacing  of  a  PhC. 

We  now  consider  the  imaging  of  an  ENGRM  slab  with  thickness  d.  The  first  surface  is  at  the 
origin  and  the  second  at  z  =  d.  A  point  source  is  placed  at  z  =  —u.  We  consider  another  point 
outside  the  ENGRM  slab  at  z  =  d  +  v.  If  the  value  of  v  satisfies  the  following  lens  equation 

u  +  v  =  ad ,  (7) 

the  phase  difference  «J>Z  =  kz(u  +  v)  +  krzd  =  Kd  is  stationary  and  an  image  will  be  formed 
without  aberration.  Note  that  for  ordinary  thin  curved  lens  with  focal  length  /,  the  lens  equation 
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Fig.  1.  Group  velocity  for  an  ENGRM  with  elliptic  EFS  and  (kr  —  icz)  •  V^cu  <  0  and 
(a)  negative  or  (b)  positive  wave  vector  refraction.  The  dashed  arrow  is  kr  —  icz.  (c)  Ray 
diagram  of  the  group  velocity  for  imaging  by  an  ENGRM  flat  lens  with  u  +  v  =  ad. 


is  h-1  +v_1  =  f~].  With  given  value  of  a  and  slab  thickness  d,  the  distance  between  the  object 
and  the  image  outside  the  slab  is  fixed,  u  +  v  +  d=  (1+  a)d.  Once  the  location  of  the  object  is 
fixed,  the  position  of  the  image  is  also  fixed,  no  matter  where  the  slab  is  placed! 

We  remark  that  a  flat  slab  can  focus  with  negative  as  well  as  with  positive  wave  vector 
refraction.  The  necessary  condition for  flat  lens  imaging  is  negative  group  refraction  .  In  general, 
there  is  not  much  meaning  of  Snell’s  law  for  wave  vectors  for  anisotropic  media.  Only  the 
refraction  law  for  the  direction  of  group  velocity  Eq.  (6)  is  meaningful  and  does  correspond  to 
ray  diagrams.  The  Veselago-Pendry  lens  is  a  special  case  of  our  flat  lens  with  a  =  1 ,  K  =  0  and 
is  the  only  case  that  the  ray  diagram  is  applicable  for  the  wave  vectors. 

The  lensmaker’s  formula  Eq.  (1)  or  (7)  only  provides  a  necessary  condition  for  imaging. 
An  additional  constraint  is  required  for  the  formation  of  a  “perfect”  image,  viz.  the  surface 
reflection  coefficient  r  should  vanish,  r  =  0,  for  far  field  and  should  diverge,  r  — >■  °°,  for  near 
field.  These  set  the  conditions  for  a  perfect  flat  lens  without  optical  axis.  The  requirement  of 
no  reflection  r  =  0  is  not  essential  for  a  slab  to  focus  far  field.  The  principal  conditions  for  the 
flat  lens  to  focus  light  are  that  Eq.  (3)  and  |  r|  -C  1  should  be  satisfied  for  all  or  at  least  a  large 
range  of  incident  angles .  The  presence  of  reflection  will  make  the  image  dim  and  may  give  rise 
to  multiple  images.  The  indefinite  indices  medium  [16]  can  have  elliptic  dispersion  relation,  but 
in  general  will  not  satisfy  the  flat  lens  imaging  conditions. 

The  field  inside  the  ENGRM  flat  lens  can  also  be  described  by  a  partial  differential  equation 
\o~2(iK+dz)2  +  d2  4-  (co/c)2] =  0.  The  medium  can  be  regarded  as  a  distorted  space  with  the 
metric  ds2  =  dx2  +  a2dz2.  Thus  there  is  an  analogy  with  gravitational  lensing  where  refraction 
occurs  due  to  the  warping  of  space  caused  by  general  relativistic  effects,  such  as  in  the  vicinity 
of  a  massive  object.  It  is  worth  noting  that  general  relativistic  effects  are  observable  due  to 
optical  refraction. 

In  the  rest  of  the  paper,  we  will  focus  on  S-polarized  electromagnetic  waves  with  the  elec¬ 
tric  field  in  the  y-direction.  In  general,  the  magnetic  permeability  p  of  the  ENGRM  should 
be  a  tensor,  and  its  relationship  to  the  other  material  parameters  a  and  k  can  be  deduced. 
Here  we  assume  that  effective  indices  e  and  p  can  be  used  to  describe  the  medium  of  in- 
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Fig.  2.  Far  field  intensity  of  a  two-pointsource  imaged  by  an  ENGRM  flat  lens  with  a  =  0.5 
and  thickness  d  =  10.  The  two  white  lines  indicate  the  surfaces  of  the  lens. 


terest.  Due  to  its  symmetry,  the  ellipsoid  can  always  be  transformed  to  its  principal  axes  [2], 
For  the  S-polarized  waves,  the  H  field  is  in  the  Jtz-plane,  only  ptx  and  nz  will  enter  our  dis¬ 
cussion.  Furthermore,  jU*  plays  a  more  active  role  than  /tz  since  the  surface  reflection  coeffi¬ 
cient  is  r  =  (jJ.xkz  —  krz)j (jixkz  +  krz).  Thus  to  achieve  perfect  imaging  for  far  field,  the  effective 
permeability  of  ENGRM  should  be 

Hx  =  K/kocos9  —  o.  (8) 

Only  in  the  case  K  =  0  that  )ix  =  —  a  is  isotropic  while  for  k  ^  0,  ju*  has  to  be  anisotropic  and 
diverges  at  kx  =  ko.  For  k  ^  0,  any  finite  nx  at  kx  =  ko  results  in  reduced  transmission.  Notice 
that  in  PhCs  including  metamaterials,  the  effective  nx  should  be  anisotropic  in  general. 

To  see  how  the  propagating  waves  are  focused,  we  consider  a  S-polarized  point  source 
— (i/4)//o^(ko|r-|-wz|)  which  is  centered  at  z  =  —  u  while  the  two  surfaces  of  the  ENGRM 
slab  are  at  z  =  0  and  z  =  d.  If  /J.x  takes  the  form  given  by  Eq.  (8),  there  will  be  no  reflection 
for  propagating  waves  and  the  transmission  coefficient  through  the  slab  is  T  =  exp (ikrzd).  The 
transmitted  far  field  is 


r n/2  dfl 

£far(r)  =  -i  exp (i Kd)  /  —  exp {ikz  [z  +  «  -  ( 1  +  a)d] }  cos (kxx) .  (9) 

Jo  2n 

Note  that  kz=ko  cos  6  and  kx  =  ko  sin  6 .  The  image  outside  the  slab  is  located  at  z  =  ( 1  +  o)d  — 
u.  One  can  see  that  except  for  a  global  phase  exp(iKtf),  the  far  field  of  the  image  is  exactly  that 
of  the  source.  The  far  field  inside  the  slab  is 


£f"r(r)  =  -/exp(iKz) 


—  exp[-/kz(az  -  «)]  cos  (kx. 


(10) 


for  0  <  z  <  d.  Another  perfect  image  is  formed  inside  the  slab  at  z  =  u/a.  The  first  surface 
acts  as  a  mirror  with  an  extra  phase  exp(j'fcz).  The  space  inside  the  ENGRM  slab  is  optically 
stretched  by  a  factor  of  a-1 .  The  flat  lens  imaging  is  shown  in  Fig.  2. 

Evanescent  waves  and  sub-wavelength  imaging  For  evanescent  waves  kx>  ko,  the  complex 
extrapolation  of  Eq.  (4)  gives 

krz  =  K+iaq  (11) 

with  q  =  (kx  —  kg)1/2.  The  transmission  coefficient  through  the  flat  lens  is  T  =  (1  — 
r1)exp(iknd)/[\  —  r2 exp(2ikrzd)\  with  r  =  [(lix-a)q  +  ik\/[(^ix  +  a)q-iK].  In  order  to  am¬ 
plify  evanescent  waves,  a  singularity  must  exist  in  T  to  compensate  to  certain  extent  the  decay 
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of  evanescent  waves  in  the  vacuum  exp[— q{u  +  v)].  In  the  case  of  single-interface  resonance 
[8, 19],  r  — ►  °°,  Texpl— q(u  +  v)\  =  exp(iKd),  all  the  evanescent  waves  will  be  amplified  only  if 

lix  =  l^  =  -o  +  iK/q.  (12) 

The  images  both  inside  and  outside  have  the  same  sub-wavelength  features  of  the  source.  If  px 
couldn’t  take  the  above  form,  sub-wavelength  imaging  is  still  possible.  At  the  so-called  overall 
resonance  condition  [19]  1  —  r2 exp(2ikrzd)  =  0,  thus  r  =  ±exp (—ikrzd),  one  gets 

px  =  p.x  =  (—o  +  iK/q)tonfffl(Gq  —  iK)d/2.  (13) 

For  k  0,  px  and  are  all  complex  with  negative  real  parts.  One  notices  that  pf  (q)  is 
flatter  than  pf (q)  as  functions  of  q.  In  terms  of  effective  indices,  the  existence  of  surface 
modes  (Oo(kx)  requires  that  91/4  be  negative  for  5-polarization.  To  amplify  evanescent  waves, 
the  c Oo(kx)  curve  must  be  very  flat  [19].  This  is  equivalent  to  say  that  the  curve  pf  (q)  should 
be  flat  so  that  the  contribution  of  evanescent  waves  are  constructive  for  certain  window  of  q. 
The  closer  px  is  to  pf  (q),  the  more  sub- wavelength  features  the  images  will  have.  Note  that 
for  large  Kd,  the  evanescent  waves  are  located  in  the  vicinity  of  the  interfaces  for  constant  px. 

Realization  in  real  materials  We  note  that  no  known  natural  material  is  found  to  have  negative 
dispersion.  Negative  group  refraction  can  be  achieved  in  periodic  or  quasi-periodic  media  with 
nonzero  K  [17].  For  most  media,  Eq.  (8)  and  (12)  are  unlikely  to  be  fulfilled,  thus  a  perfect  flat 
lens  is  unattainable.  Realization  using  real  materials  is  discussed  next. 

Ordinary  material  may  not  have  the  EFS  described  by  Eq.  (4)  for  all  kx.  For  an  uniaxial  crys¬ 
tal,  since  the  group  refraction  is  positive,  thus  krz  =  K+\o\kz  with  k  =  0.  There  will  be  no  focus 
though  negative  refraction  for  certain  incident  angles  can  be  easily  achieved  [18].  Instead  a  vir¬ 
tual  image  will  be  formed  satisfying  the  same  lens  equation.  However,  a  PhC  could  have  both 
an  elliptic  EFS  and  negative  group  refraction  for  certain  windows  of  kx  at  certain  frequency. 
Here  we  use  some  general  features  of  PhC  to  explain  the  imaging  mechanism.  To  this  end,  we 
consider  the  first  band  of  a  square  lattice  PhC  of  unit  lattice  spacing.  For  small  kx  along  the  TM 
direction,  the  dispersion  can  be  approximated  elliptically  as  krz  —  K  —  Ookz  with  appropriate 
constants  k  and  ao  (see  Fig.  3).  At  too  =  (Ou,  the  upper  limit  of  all-angle  negative  refraction 
in  the  first  band  [12],  k  =  \fln  —  &o(l  —  (To).  To  enhance  the  transmission,  the  slab  thickness 
can  be  selected  to  satisfy  the  Fabry-Perot  resonance  condition  (k  —  Oo(o/c)d  =  0  (mod  2k). 
The  slab  is  on  resonance  only  for  kx  — >  0.  For  large  kx,  krz  ~  \[2n  —  a0  l^2kz,  large  transmis¬ 
sion  would  require  px  to  diverge.  Thus  the  image  location  is  bounded  Ob  <  (u  4-  v)/d  <  <70 
and  the  Fourier  components  of  the  object  with  kx  — +  ko  will  be  partially  lost.  This  leads  to  a 
so-called  self-collimation  effect  for  small  (To-  Far  field  images  both  inside  and  outside  the  slab 
with  reasonably  good  quality  will  be  formed.  However,  the  image  inside  the  PhC  is  stretched 
and  modulated  due  to  the  partial  transmission  besides  the  Bloch  modulation.  Even  with  the  re¬ 
moval  of  Bloch  modulation,  the  image  inside  PhC  will  be  visible  only  for  large  thickness.  For 
evanescent  waves  at  w  =  tou  with  kx  — >  ko,  one  has  ^px  =  —  cr0  and  3px  — >  °°  while  for 
kx  >•  ko,  px(kx)  ~  Hx(kx)  -+  —  (Tjj"1.  To  have  substantial  sub-wavelength  feature  of  the  image, 
the  real  part  of  fix  must  take  large  negative  value  for  small  Oo.  An  example  is  shown  in  Fig.  3, 
which  is  very  similar  to  Fig.  5  in  Ref.  [12].  The  image  is  confined  in  the  vicinity  of  the  second 
surface  of  the  PhC  slab  due  to  a  small  value  of  Ob. 

In  this  paper,  we  have  discussed  the  group  refraction  of  an  anisotropic  medium  characterized 
by  two  materials  parameters  a  and  k.  When  a  is  a  positive  constant,  the  anisotropic  medium 
is  ENGRM  and  a  flat  slab  of  such  medium  can  be  used  as  a  focal  lens  without  optical  axis  and 
leading  to  images  free  of  aberration.  Our  theory  of  flat  lens  imaging  is  a  generalization  of  the 
Vesalago-Pendry  perfect  lens  and  beyond  Silin’s  formula.  The  theory  is  valid  for  real  and  virtual 
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Fig.  3.  (a)  EFS  of  the  TE  modes  of  a  square  lattice  PhC  calculated  using  plane  wave  ex¬ 
pansion.  (b)  Far-field  Hz  of  a  point  source  across  an  ENGRM  slab  with  dispersion  relation 
kn  =  K—okz  with  K  =  3.2465  and  a  =  0.08  which  approximates  the  actual  EFS  of  the 
PhC  (dashed  line  in  the  left  panel).  The  permittivity  ex  =  2.2  is  used.  Note  the  modulated 
field  inside  the  homogeneous  ENGRM.  Due  to  impedance  mismatch,  the  images  are  not 
“perfect”  unlike  Fig.  2. 


images.  This  flat  lens  theory  leads  to  a  clear  understanding  of  negative  group  refraction  and 
flat  lens  imaging  in  electromagnetism  [13, 12],  acoustics  [20],  and  electron  waves.  Numerical 
simulations  of  homogeneous  ENGRM  and  of  PhCs  were  carried  out  supporting  our  theory. 

The  lack  of  optical  axis  for  the  flat  lens  has  very  broad  applications  and  confers  important 
advantages  in  optics.  Clearly  there  are  no  aperture  restrictions.  There  are  two  characteristics 
that  might  be  viewed  as  limitations.  The  lensmaker’s  formula  clearly  dictates  upper  bounds  to 
u,v  <  ad  for  real  image  formation.  The  other  limitation  is  that  magnification  is  always  unity. 
For  any  medium  the  lens  properties  a,  jc,  and  also  the  working  frequency  co o  can  be  obtained 
by  inspection  of  EFS.  To  have  high  quality  image,  a  should  be  constant  and  close  to  unity. 
The  present  theory  can  be  used  to  design  tailor-made  flat  lens.  Extension  of  our  theory  to  three- 
dimensions  is  straightforward. 

Work  supported  by  the  National  Science  Foundation  and  the  Air  Force  Research  Laborato¬ 
ries,  Hanscom,  MA. 
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ABSTRACT 

Negative  refraction  and  left-handed  electromagnetism  in  a  photonic  crystal  are  demonstrated  in  waveguide  and 
free  space  experiments  at  microwave  frequencies.  Precision  control  to  achieve  tailor-made  refractive  indices  has  been 
achieved.  The  negative  refraction  in  these  photonic  crystals  is  shown  to  lead  to  imaging  by  a  flat  lens.  We  have  also 
developed  a  generalized  theory  of  flat  lens  imaging.  These  results  promise  potential  applications  in  a  variety  of  optical 
and  microwave  systems  for  communications  and  imaging. 


1.  Introduction 

Negative  refraction  (NR)  is  a  new  phenomenon  allowed  by  Maxwell’s  equations  that  has  led  to  novel  effects 
regarding  the  propagation  and  control  of  electromagnetic  waves.  In  a  negatively  refracting  material,  electromagnetic 
waves  obey  a  left-handed  relationship  between  the  vectors  E,  H  and  k,  and  hence  the  media  are  also  referred  to  as  left- 
handed  metamaterials  (LHM).  NR  has  been  recently  demonstrated  in  two  types  of  LHM  -  composite  metamaterials 
made  of  split  ring  resonators  and  wire  strips,  and  in  photonic  crystals  (PC) 1V,V. 

Negative  refraction  in  LHM  leads  to  new  optical  and  electromagnetic  components  that  open  the  door  for  new 
applications  from  microwave  to  optical  frequencies.  An  interesting  application  of  negative  refraction  is  superlensingvl 
effect  by  a  flat  lens  with  no  curved  surfaces  that  can  potentially  overcome  the  diffraction  limit  imposed  by  conventional 
lenses.  We  have  recently  fabricated  a  flat  lens  recently  using  a  photonic  crystal  structure  and  demonstrated  focussing 
without  optical  axis™.  The  flat  lens  is  best  suited  when  the  source  is  close  to  the  surface  of  the  lens.  For  those 
applications  where  imaging  of  far  objects  is  important,  we  have  recently  demonstrated  focusing  by  a  plano-concave  lens 
using  negative  refraction,  thus  enabling  several  applications  of  far  field  imaging™1. 

In  this  paper  we  describe  the  experimental  setups  and  results  that  demonstrate  negative  refraction  in  photonic 
crystals.  The  experimental  results  are  in  excellent  agreement  with  band  structure  calculations.  A  generalized  theory  of 
imaging  by  a  flat  lens  is  briefly  described,  accompanied  by  numerical  confirmation  of  imaging  by  a  lens  composed  of 
periodic  scatterers.  These  results  confirm  the  experimental  results  reported  earlier  in  real  photonic  crystal  flat  lens 
experiments. 

2.  Experimental  setups 

We  have  carried  out  two  types  of  experiments  that  demonstrate  negative  refraction  in  PC,  parallel  plate 
waveguide'x,x  and  free  space  experiments  inside  an  anechoic  chamber111.  In  both  cases  microwave  refraction 
measurements  were  carried  out  on  wedge  prisms  fabricated  from  PC  structures,  and  negative  refraction  was  clearly 
observed. 

2.1.  Parallel  Plate  Waveguide  experiments 

The  microwave  photonic  crystal  consists  of  an  array  of  cylindrical  copper  rods  with  a  radius  0.63  cm  arranged 
in  a  triangular  lattice.  The  ratio  of  the  radius  r  to  lattice  constant  a  is  kept  atr/a  =  0.2  .  Refraction  measurements  are 
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Figure  1  (Left)  Free  space  experimental  setup  for  microwave  refraction  measurements.  (Right) 
Directions  of  the  wavevectors  and  the  group  velocity  in  the  refraction  experiments  for  negative  and 
positive  refraction. 


carried  out  in  a  metallic  parallel  plate  waveguide.  The  excitation  inside  this  waveguide  is  transverse  magnetic  (TM), 
meaning  that  the  electric  field  Ey ,  is  parallel  to  the  axis  of  the  dielectric  rods.  The  collimated  microwave  beam 

propagates  a  distance  of  ~lm  in  the  waveguide,  and  then  is  incident  normally  to  the  PC  wedge  (see  Fig.  1). 

The  other  face  of  the  prism,  which  is  the  surface  of  refraction,  is  at  60°  to  the  incident  beam.  A  dipole  antenna 
is  mounted  on  a  goniometer  that  runs  along  the  semicircular  outer  edge  of  the  parallel  plate  waveguide  to  detect  the 
refracted  beam.  The  refraction  is  considered  positive  (negative)  if  the  emerging  beam  is  detected  at  positive  (negative) 
angles  as  shown  in  Fig  1. 

The  calibration  of  the  refraction  experiment  set  up  is  validated  by  the  refraction  experiment  performed  on  a 
polystyrene  prism  of  the  similar  dimensions  to  the  PC  prism  that  corresponds  to  a  positive  refractive  index  n  =  1.58  as 
expected  by  theory. 


2.2.  Free  Space  Experiments 

In  addition  to  the  parallel  plate  experiments  described  above,  we  also  carried  out  free  space  (FS)  experiments 
using  metallic  PC  prism  suspended  in  air.  The  principal  advantage  of  these  free  space  experiments  is  the  ability  to 
separate  out  beams  over  large  transmission  distances  of  3m.  Furthermore,  the  polarization  can  be  varied  and  refraction 

experiments  were  carried  out  for  both  for  both  TM  (  E  \\  to  the  rod  axis)  and  TE  (E  _L  to  the  rod  axis)  mode 
propagation.  Excellent  agreement  with  band  structure  calculations  was  observed  for  both  polarizations,  showing  that  a 
PC  can  exhibit  negative  refraction  with  tailor  made  refractive  indices  in  a  large  frequency  range. 

The  microwave  photonic  crystal  consists  of  an  array  of  cylindrical  copper  tubes  arranged  on  the  same 
periodical  triangular  lattice  as  described  in  the  above  parallel  plate  waveguide  experiment.  Refraction  experiments  were 
performed  in  an  anechoic  chamber  of  dimensions  5  x  8  x  4  m3  to  prevent  reflections  from  the  walls.  A  square  X-band 
horn  placed  at  3  m  from  the  PC  acts  as  a  plane  wave  source  (Fig.l).  A  piece  of  microwave  absorber  with  a  6x6inch2 
aperture  in  front  of  the  PC  was  used  to  narrow  the  incident  beam.  On  the  far  side  another  square  hom  was  attached  to  a 
goniometer  which  swings  around  in  two-degree  steps  to  receive  the  emerging  beam.  Refraction  is  considered  positive 
(negative)  if  the  emerging  signal  is  received  to  the  right  (left)  of  the  normal  to  the  surface  of  refraction  of  the  PC. 

Measurements  were  carried  out  with  the  incident  wave  vector  ki  along  r  ->  M  (0,1)  direction  of  the  first  Brillion  zone 
of  the  PC  and  in  both  TM  and  TE  modes.  The  angle  of  incidence  0  =  60°  for  r  M  is  chosen  in  order  to  minimize 
surface  periodicity  along  the  surface  of  refraction,  thus  eliminating  higher  order  Bragg  waves. 


3.  Results 

A  3D  plot  of  the  transmitted  intensity  measured  at  different  angles  in  the  frequency  window  6-12GHz  for  the  TM  mode 
propagation  is  shown  if  Fig.  2  b2.  Also  shown  are  polar  plots  of  the  transmitted  intensity  at  selected  frequencies  6.68, 
9.315  and  10.77  GHz  showing  the  emerging  beams. 
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As  can  be  seen  from  the  figure  between  6  and  7.1  GHz  the  signal  emerges  on  the  positive  side  of  the  normal  to  the 
surface  corresponding  to  positive  refraction.  No  transmission  is  observed  between  7.1  and  8.3  GHz.  Above  8.3  GHz  and 
1 1  GHz  two  signals  are  observed  on  the  positive  and  negative  sides  of  the  normal  corresponding  to  1st  and  0th  order 
Bragg  wave  refraction.  The  negatively  refracted  signal  is  strongest  from  10.14-11  GHz  and  positively  refracted  signal 
around  8.6-9.6  GHz  and  in  the  frequency  window  9.6-10.14  GHz  both  signals  have  the  same  strength.  So  with  the 
increase  in  frequency  the  positive  signal  gets  weaker  while  negative  signal  gets  stronger.  The  experimental  refractive 
index  n  is  obtained  from  applying  Snell’s  law  n  =  sin(0r)/sin(0;)to  each  emerging  beam.  The  validity  of  Snell’s  law 
has  been  established  earlier  in  metallic  PCs  (ix). 


f  =  6.68  GHz 

0  al) 


f  =9.315  GHz  a2) 

0 


Fig.  2  Figures  2  al,  a2,  a3  show  polar  plots  of  wave  directivity  for  representative 
frequencies  of  6.68  GHz,  9.315  GHz  and  10. 77GHz.  Figure  3bl  shows  the  band  structure  of 
the  triangular  lattice  of  metallic  PC  with  packing  density  r/a=0.2  for  TM  Mode  and  Figure 
2b2  shows  the  3D  experimental  map  of  the  refracted field  versus  frequency  and  angle  of 
refraction. 


18 


3 


Refractive  angle 


Fig.  3  3a  shows  a  3D  experimental  plot  of  the  magnitude  of  the  field  versus  frequency  and  angle  of  refraction 
and  Fig.  3  b  shows  the  band  structure  calculation  for  triangular  lattice  PC  with  r/a=0.2  in  TE  mode 


We  have  also  carried  out  measurements  of  refraction  for  TE  mode  propagation.  The  band  structure  for  TE 
mode  is  shown  in  the  figure  3.  Here  negative  refraction  is  observed  between  6.4  -  9.8  GHz  and  positive  refraction 
between  9.8  -12  GHz.  It  is  important  to  note  that  negative  refraction  is  possible  for  both  TM  and  TE  modes;  such  a 
freedom  in  the  choice  of  modes  provides  a  crucial  advantage  of  using  the  metallic  PC  over  the  split  ring  wire  array 
metamaterial. 


3.1.  Band  Structure  Calculations 

The  physical  principle  behind  the  present  results  can  be  understood  from  the  band  structure  of  the  metallic  PC. 
We  have  calculated  the  band  structure  of  the  triangular  lattice  PC  employing  a  finite  difference  method.  For  the  TM 
modes,  the  electric  field  can  be  identified  with  the  wave  function  of  a  2D  periodic  8-shell  potential  barrier,  thus  a  plane 
wave  expansion  method  can  be  used  to  calculate  the  band  structure.  The  2D  band  structures  of  both  TM  and  TE  are 

shown  respectively  in  Figs.  2  (bl)  and  3(b).  For  a  plane  wave  with  incident  wave  vector  k-  and  frequency  CO  incident 

normally  on  an  air-PC  interface,  the  wave  vector  kj  inside  the  PC  is  parallel  or  anti-parallel  to  kt  as  determined  by  the 

band  structure.  For  a  general  case  the  phase  velocity  vp=(c/\np\)kf  with  kf=kft |fcy|  and  the  group  velocity  vg  =  V-tu 

can  have  any  angle  between  them.  Thus  distinction  must  be  made  between  negative  wave  vector  refraction  and  negative 
group  refraction. 

An  effective  refractive  index  can  be  defined  np  =  sgn(vK  kf  1  /  co  and  calculated  from  the  band  structure. 

The  sign  of  np  is  determined  from  the  behavior  of  the  EFS.  Negative  wave  vector  refraction  will  happen  if 

vg-kf  <  0  while  negative  group  refraction  will  be  present  only  if  vg-(kf  -k)<  0  with  k  is  the  center  of  the  equi- 
frequency  surface  (EFS)  ellipsoid.  Note  that  the  direction  of  the  group  velocity  vg  in  an  infinite  PC  coincides  with  that 

of  the  energy  flow.  If  only  the  zero-th  order  Bragg  wave  refraction  is  present,  conservation  of  the  kj-  component  along 
the  surface  of  refraction  would  result  in  positive  or  negative  wave  vector  refraction,  depending  on  whether  Ay  is  parallel 
or  anti-parallel  to  group  velocity. 

If  higher  order  Bragg  wave  refraction  is  present,  the  emerging  beam  can  be  written  as  'P,  =  a0e,k,°'r  +  a{e,k'yr 
where  kt o  and  An  represent  the  refracted  wave  vectors  corresponding  to  the  0th  and  lsl  order  Bragg  wave-vectors  of  the 
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field  inside  the  PC.  2nd  and  higher  order  Bragg  wave  refraction  is  ignored.  Let  the  component  of  kt  along  the  normal  to 
the  surface  of  refraction  of  the  prism  be  ktl .  One  has  k,0±  =  ^m2 /c2 -kj  sin2  0  and*(]±  =  ^a2/c2 -(2x/a-kfsin0)2 
where  Q=n/3  is  the  incident  angle  at  the  second  surface  of  the  prism  and  a  is  the  lattice  periodicity.  The  parallel 
component  is  =  kj  sin  0 .  The  refracted  angle  can  be  obtained  from  0r  =  tan-1  (A^  / k,L)  for  each  beam. 

From  the  band  structure  and  the  EFS  for  TM  mode  of  propagation  negative  refraction  is  predicted  for  the 
second  and  third  band  regions,  with  positive  refraction  in  the  1st  band.  In  the  1st  band  between  6-7.1  GHz  the  EFS  move 
outward  with  increasing  frequency,  so  that  vg  -  kf  >  0-  In  the  second  band  between  8.3-11  GHz,  the  EFS  move  inward 

with  increasing  frequency,  consistent  with  np<  0  corresponding  to  vg  ■  Ay  <  0  (  vg  anti-parallel  to  kj  ).  The  band  gap  is 

in  the  frequency  range  7. 1-8.3  GHz  between  the  1st  and  2nd  pass-bands  and  from  1 1-1 1.2  GHz  between  the  second  and 
third  bands. 

3.2.  Discussion 
TM  MODE 

The  agreement  between  experimental  results  of  refraction  and  those  deduced  from  band  structure  calculation  is 
remarkable.  In  Fig  3(c)  the  representative  experimental  results  show  that  negative  and  positive  refraction  are  present  in 
different  frequency  ranges. 

I.  In  the  first  band  between  6-7.1  GHz  the  refraction  is  positive,  consistent  with  np>  0  from  the  band  structure 

and  corresponding  to  RHE  vg  •  Ay  >  0 .  The  representative  plot  at  f  =  6.68  GHz,  Fig.  2(al)  shows  that  the  signal 
is  largely  in  the  positive  direction. 

II.  In  the  second  band  between  8.3  -  1 1  GHz,  the  refraction  is  negative  for  the  zero’th  order  Bragg  wave  refraction, 
consistent  with  np  <  0  from  the  band  structure  and  corresponding  to  LHE  vg  -  kf  <  0  .There  is  observed  also  a 

first  order  Bragg  wave  that  is  positively  refracted  in  the  same  frequency  window.  This  branching  of  the  beams 
observed  in  metallic  rods  in  air  design  can  be  eliminated  in  a  dielectric  PC  design,  where  only  a  single  beam  can 
be  observed.  Two  illustrative  plots  corresponding  to  9.315  and  10.77  GHz  show  that  the  signal  is  strongly  peaked 
in  the  negative  direction  and  becomes  stronger  with  the  increase  in  frequency. 

III.  A  band  gap  is  clearly  observed  in  the  experiment,  which  matches  with  the  band  structure,  in  the  region  between 
7. 1  -8.3  GHz.  As  can  be  seen  from  the  3D  plot  figure  the  transmission  through  the  PC  is  near  zero. 


TE  MODE 

It  is  important  to  note  that  negative  refraction  is  possible  for  TE  mode  as  well.  Such  a  freedom  in  the  choice  of  modes  is 
a  crucial  advantage  of  using  the  metallic  PC  over  the  split  ring  wire  array  metamaterial. 

I.  In  the  first  band  between  6.4  -  9.8  GHz  the  zero’th  order  Bragg  wave  refraction  is  negative,  consistent  with 
np<  0  and  vg  -  kf  <  0  from  the  band  structure  and  corresponding  to  LHE.  The  bandwidth  of  this  negative  band 
is  39%  which  is  quite  a  remarkable  improvement. 
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QM  square  lattice,  r/a=0.2,  y=-4.5  (5-potential) 


Real,  to/2jr=0.14,  u =0.396 45 


Flat  lens  image  of  a  periodic 
delta-shell  potential 


The  intensity  of  the  transmitted 
wave  function. 


Fig.  4  a)  Band  structure  of  a  periodic  delta- 
shell  potential  of  a  square  lattice.  The  straight 
lines  are  light  lines  b)  Flat  lens  image  of  a 
periodic  delta-shell  potential  c)  The  intensity 
of  the  transmitted  wave  junction. 


II.  In  the  second  band  between  9.8-12  GHz,  the  refraction  is  largely  positive,  consistent  with  np>  0  from  the 
band  structure  and  corresponding  to  RHE  vg-k/>  0 . 

III.  The  absence  of  a  band  gap  from  6-1 2GHz  along  GM  Brillion  zone  propagation  direction  as  predicted  by  the  band 
structure  calculation  is  in  complete  agreement  with  experimental  results  shown  in  Fig  3a. 


3.3.  Imaging  by  a  Flat  Lens 

Negative  refraction  using  LHM  leads  to  entirely  new  and  fascinating  concepts  for  image  formation.  One  of  the  most 
striking  consequences  of  NR  is  the  ability  of  flat  surfaces  to  focus  a  divergent  beam  and  hence  leads  to  imaging  of  an 
object  placed  near  the  lens.  A  negative  index  flat  lens  reconstructs  not  only  the  far  field  propagating  waves,  but  can  also 
reconstruct  evanescent  waves. 
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Image  formation  is  wave  reconstruction  of  both  far  and  near  fields.  We  have  developed  a  new  approach  to  imaging  by  a 
flat  lens  made  of  a  proposed  artificial  metamaterial  [xii].  The  key  aspect  of  a  flat  lens  using  LHM  is  the  absence  of  an 
optical  axis.  In  contrast,  current  flat  lenses  like  the  Fresnel  lens,  which  use  gradient-index  refractive  material  (GRIN), 
and  all  have  an  optical  axis.  The  conditions  for  image  formation  are  determined  by  the  requirement  that  the  phase  shift 
<D  =  Jit  ■  dl  from  object  to  image  is  independent  of  wave  vector  for  all  directions.  Flat  lens  imaging  can  be  achieved  by 

fabricating  a  material  which  has  appropriate  dispersion  relations  m(k) .  The  intrinsic  parameters  that  define  the  new 

material  are  o;  which  is  determined  by  the  dispersion  characteristics  of  the  material,  and  k,  which  is  a  measure  of  the 
phase  shift  exp(/cr«/)  by  the  flat  slab.  This  is  a  generalization  of  Pendry’s  flat  lens  with  elliptic  dispersion 

k\  +  cr~2(kz  -  tcj1  -of  /  c2  and  negative  group  refraction  vg  ■  (kf  -  k)  <  0 .  The  lens  equation  for  a  slab  of  thickness  d  is 
determined  to  be  u+v=ad.  Here  u  and  v  are  the  distances  of  the  object  and  the  image  from  the  nearest  surfaces.  Unlike 
GRIN,  the  effective  refractive  index  is  angle  dependent  n{9) ,  not  spatial  dependent  n(x) .  The  lens  constant  k  determines 
whether  the  lens  has  positive  or  negative  refractive  index  np  and  gives  rise  to  a  phase  shift  of  the  image  exp(z'er/a/) . 
Simulations  of  image  formation  by  such  a  material  are  shown  in  Fig.7.  The  space  inside  the  slab  appears  to  be  stretched 
by  a  factor  of  a.  “Perfect”  images  are  achieved  when  the  material  is  suitably  matched  to  its  environment  as  shown  in  the 
figure.  We  have  shown  that  this  material  can  be  realized  approximately  using  photonic  crystals  to  achieve  the  desired 
a>(k)  .This  preliminary  theoretical  investigation  of  a  special  material  for  flat  lens  imaging  is  entirely  novel  and  can  lead 
to  new  designs  for  flat  lenses.  The  theoretical  principles  that  we  have  obtained  will  be  used  to  design  novel  flat  lenses 
made  of  metamaterials  and  photonic  crystals.  The  main  goal  is  to  design  a  material  with  suitable  characteristics  leading 
to  minimal  reflection  and  “perfect”  imaging. 

In  our  flat  lens  experiment  the  image  is  imperfect  and  reveals  the  presence  of  satellites  -  this  is  due  to  the  fact  that  the 
propagation  is  along  the  so-called  T  -»  X  for  a  square  lattice  in  the  2nd  band  where  the  inside  the  crystal  the  1st  Bragg 
wave-vector  is  involved.  A  better  image  should  be  achieved  with  propagation  along  the  T  ->  M  in  the  first  band  where 
also  the  negative  group  refraction  condition  can  be  satisfied.  Under  this  condition  a  single  image  should  be  possible. 
More  generally  the  best  conditions  of  imaging  are  achieved  when  the  equifrequency  surfaces  are  circular  or  nearly  so, 
and  consistent  with  negative  refraction. 

NUMERICAL  SIMULATIONS  OF  FLAT  LENS  IMAGING 

To  understand  flat  lens  imaging,  we  consider  a  model  system  which  is  a  periodic  8-shell  potential 

V(r)=  /y'.jSQr  —  R j  \  -r0) .  For  y— >oo,  the  wave  function  must  vanish  on  the  disk  boundary.  The  wave  function  is 

equivalent  to  the  electric  field  of  the  TM  modes  of  a  metallic  PC.  We  have  used  this  model  to  simulate  negative 
refraction  in  metallic  PC.  For  attractive  potential  y  <  0,  and  within  appropriate  energy  window,  negative  refraction  and 
flat  lens  imaging  of  waves  can  be  achieved.  Below  is  a  band  structure  of  a  square  lattice  for  y=-4.5  and  rfa=0.2  with 
lattice  spacing  a.  It  is  predicted  that  above  /ta/2^=0.14,  all-angle-negative-refraction  can  be  achieved  along  the  TM 
direction. 

The  flat  lens  imaging  of  a  point  source  -±H0(kr)  is  shown  below  at  k  a/2n=  0.14  with  a  distance  u  —  0.4a  away  from 

the  first  surface.  A  clear  image  is  seen  on  the  other  side  of  the  slab  at  v=4.70a.  This  gives  an  effective  aefl=1.03  which  is 
consistent  with  the  EFS.  When  the  source  is  moved  away  from  the  first  surface  of  the  slab,  the  image  will  move  toward 
the  second  surface.  When  u  >  ad,  virtual  image  will  be  seen  through  ray-tracing. 

The  mechanism  of  sub-wavelength  imaging  is  due  to  the  bounded  states  of  the  8-shell  potential.  The  surface  waves  are 
not  present  if  no  8-shell  is  broken. 

For  real  photonic  crystals  or  the  above  8-shell  potential  model,  the  effective  refractive  index  a  is  angle  dependent.  At 
certain  frequency  and  with  appropriate  design,  the  EFS  can  be  elliptic  with  negative  group  refraction,  the  image  can 
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have  high  quality  with  small  aberration.  Otherwise,  the  image  will  be  blurred  due  to  aberration.  Our  flat  lens  theory  can 
be  used  to  design  high  quality  flat  lens. 

4.  Conclusion 

The  phenomenon  of  negative  refraction  in  periodic  media  embodied  by  photonic  crystals  is  well  established  by  the 
experiments  described  above.  The  imaging  and  focusing  experiments  also  show  that  entirely  novel  concepts  are  possible 
due  to  NR.  Optical  components  utilizing  NR  are  expected  to  have  several  advantages  such  as  enabling  sub-wavelength 
resolution,  flat  lenses  without  optical  axis  thus  avoiding  optical  alignment,  and  improved  performance  such  as  due  to 
reduced  aberration  for  the  same  structural  dimensions.  In  addition  light  weight  and  compact  structures  offer  additional 
advantages  in  a  variety  of  applications. 

We  thank  Beverly  Turchenetz  and  Richard  Wing  for  invaluable  contributions.  This  work  was  supported  by  the 
National  Science  Foundation  and  the  Air  Force  Research  Laboratories,  Hanscom  AFB  (contract  #F33615-01-1-1007). 
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We  demonstrate  focusing  of  a  plane  microwave  by  a  planoconcave  lens  fabricated  from  a  photonic 
crystal  having  a  negative  refractive  index  and  left-handed  electromagnetic  properties.  An  inverse 
experiment,  in  which  a  plane  wave  is  produced  from  a  source  placed  at  the  focal  point  of  the  lens, 
is  also  reported.  A  frequency-dependent  negative  refractive  index,  n(to)  <  0  is  obtained  for  the  lens 
from  the  experimental  data  which  match  well  with  that  determined  from  band  structure 
calculations.  ©  2005  American  Institute  of  Physics.  [DOI:  10.1063/1.1927712] 


Negative  refraction  in  left-handed  materials  (LHM)'“6 
has  triggered  intense  interest  in  designing  microwave  and 
optical  elements,  a  flat  lens  being  one  of  them.7  In  the  case  of 
the  flat  lens,  the  waves  entering  from  the  source  refract  nega¬ 
tively  on  both  interfaces  and  meet  constructively  on  the  far 
side  of  it.  Thus,  a  flat  lens  applies  phase  correction  to  the 
propagating  waves — similar  to  a  conventional  lens  made  of  a 
naturally  available  material  and  having  a  positive  index  of 
refraction.  However,  it  operates  only  when  the  source  is 
close  to  the  lens.7,8  For  a  majority  of  applications  of  lenses  in 
optics,  astronomical  telescopes,  commercial,  and  defense  mi¬ 
crowave  communications,  far-field  imaging  is  required. 
Negative  refraction  allows  focusing  of  a  far-field  radiation  by 
concave  rather  than  convex  surfaces,9  with  the  advantage  of 
reduced  aberration  for  the  same  radius  of  curvature  and 
changes  many  commonly  accepted  aspects  of  conventional 
optical  systems.  Of  the  two  classes  of  LHM  currently  being 
investigated,  focusing  using  a  planoconcave  lens  made  of  a 
left-handed  metamaterial  (MM)  fabricated  by  interleaving  ar¬ 
rays  of  wire  strips  and  split  ring  resonators  was  demonstrated 
experimentally  by  Parazzoli  et  al.10 

In  this  letter,  we  demonstrate  that  a  real  image  of  a  far- 
field  radiation  can  be  produced  using  a  left-handed  photonic 
crystal  (PhC)  lens.  We  also  report  an  inverse  experiment  in 
which  the  lens  produces  plane  waves  from  a  point  source 
placed  at  the  focal  length.  The  frequency-dependent  refrac¬ 
tive  index,  n{a>),  determined  from  the  experimental  data,  is 
in  complete  agreement  with  that  predicted  by  the  theory  us¬ 
ing  band  structure  calculations.  The  results  confirm  that  far- 
field  focusing  is  realizable  and  opens  the  door  for  several 
applications  of  the  LHM  in  the  far-field  region. 

Microwave  focusing  measurements  are  carried  out  using 
three  planoconcave  lenses  made  of  a  dielectric  PhC.  The 
radii  of  curvature  of  the  lenses  are  13.5,  17.5,  and  22  cm. 
The  two-dimensional  PhC  consists  of  a  periodic  array  of 
alumina  rods  in  air,  arranged  on  a  square  lattice  and  having 
dielectric  constant,  e=8.9.  The  ratio  of  the  radius  of  the  alu¬ 
mina  rods  to  the  lattice  spacing  is  r/a=0.175.  Microwave 
measurements  are  carried  out  in  a  parallel  plate  waveguide. 
An  Y-band  waveguide  kept  at  a  distance  of  150  cm  from  the 
flat  surface  of  the  lens  acts  as  a  microwave  source.  The  emit¬ 
ted  wave  travels  through  the  parallel-plate  waveguide  and 
the  eventual  plane  wave  is  made  to  incident  on  the  flat  sur¬ 


face  of  the  lens.  The  propagation  of  the  wave  inside  the  PhC 
is  along  I"Y  direction  of  the  Brillouin  zone.  Field  maps  of  the 
incoming  plane  wave  and  the  emerging  radiation,  on  the  far 
side,  are  captured  using  a  monopole  sensor  on  a  ground 
plane.  The  sensor  is  hooked  up  to  an  automated  X-Y  transla¬ 
tional  stage  which  scans  for  the  electric-field  component  of 
the  microwaves  in  the  region  of  interest  as  shown  in  Fig.  1 . 
An  HP-85 10C  network  analyzer  is  used  for  measuring  the 
transmission  characteristics. 

Figure  2  shows  a  sharp  focus  point  achieved  at  10.1  cm 
using  the  planoconcave  PhC  lens  of  radius  of  curvature  13.5 
cm.  From  left  to  right,  the  incoming  plane  wave,  a  real  pic¬ 
ture  of  the  PhC  lens,  and  the  emerging  mapped  field  are 
shown.  A  clear  focusing  point  is  observed  in  the  frequency 
range  of  9.265-9.490  GHz.  Note  that  the  direction  of  the 
energy  flow  changes  only  at  the  second  interface  of  the 
planoconcave  lens.  To  validate  that  the  focusing  is  due  to 
negative  refraction,  an  inverse  experiment  is  carried  out,  in 
which  a  point  source  is  kept  at  the  observed  focal  point  of 
the  lens.  As  can  be  seen  from  Fig.  3,  a  circular  wave  front 
from  the  point  source  after  passing  through  the  lens  emerges 
as  a  plane  wave.  These  two  results  validate  the  behavior  of 
the  left-handed  planoconcave  lens.  A  similar  lens  made  of  a 
naturally  available  material,  on  the  contrary,  would  not  be 
able  to  produce  a  real  focus  point  as  the  present  PhC  lens 
does,  but  rather  produce  a  diverging  beam. 

The  refractive  index  of  the  lens  is  determined  using  the 
lens  equation  n=l  -Rtf,  where  R  is  the  radius  of  curvature 


Scanned  area 
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FIG.  2.  (Color)  Focusing  by  a  planoconcave  PhC  lens  having  radius  of 
curvature  13.5  cm.  The  focus  point  observed  at  9.31  GHz  is  10.1  cm  from 
the  concave  lens  surface.  A  photograph  of  the  PhC  is  superimposed  on  two 
matlab  surface  plots  to  obtain  the  final  figure.  Dark  strip  in  the  center  is  a 
schematic  representation  of  the  area  between  the  lens  and  the  incoming 
wave.  On  the  left  side,  field  map  of  the  incoming  plane  wave  is  shown  (real 
part  of  transmission  coefficient)  and  on  the  right  side,  intensity  of  the  focus 
point.  Scale:  On  the  left,  the  real  part  of  S21  varies  from  -0.025  to  0.025,  on 
the  right  side  intensity  from  0  to  1.6  X  10~3.  Dimensions  of  the  figure  are 
49  X  34  cm2.  PhC  lattice  spacing  is  1.8  cm  and  the  packing  density  of  the 
square  lattice  is  determined  from  the  ratio  r/a=0.175. 

and  /  is  the  focal  length.  Using  this  description,  we  get 
n=-0.4  at  9.25  GHz.  Note  that  real  focusing  by  a  planocon¬ 
vex  lens  is  achieved  with  n>  1 ,  R< 0,  while  for  the  plano¬ 
concave  lens  with  R>0,  n<  1.  In  Fig.  4,  experimentally  de¬ 
termined  n(co)  for  all  three  lenses  is  shown.  It  can  be  seen 
that  for  R=  13.5  cm  a  sharp  focus  is  achieved  in  the  fre¬ 
quency  range  from  9.25  GHz  to  9.5  GHz,  for  R=  17.5  cm  in 
the  range  from  8.5  GHz  to  9.2  GHz,  and  for  R= 22  cm,  from 
7.8  GHz  to  8.1  GHz. 

The  nature  of  the  left-handed  electromagnetism  and  fo¬ 
cusing  can  be  understood  from  the  dispersion  characteristics 
of  the  PhC.  Figure  5  shows  band  structure  calculated  for  the 
PhC  using  a  plane  wave  expansion  method.  From  the  band 
structure,  it  can  be  deduced  that  in  the  second  band  for 
propagation  along  FA-  direction,  the  wave  vector  k  is  in  op¬ 
posite  direction  to  group  velocity,  vg:vg-k<0,u'5  resulting 
in  negative  refraction  in  the  second  band  and  correspond¬ 
ingly  negative  refractive  indices.  From  the  band  structure  the 
theoretical  refractive  index  n=ck/(o  is  determined.  The  solid 
line  in  Fig.  4  shows  n{a>)  obtained  at  various  frequencies.  As 
the  frequency  increases  from  7.8  to  9.8  GHz,  the  theoretical 
n  increases  from  -1.1  to  0.  Note  that  both  the  experimentally 
obtained  refractive  indices  and  the  theoretically  calculated 
values  are  in  complete  agreement. 

In  Fig.  4,  it  can  be  observed  that  n(a>)  determined  for 
each  radius  of  curvature  fall  in  certain  frequency  windows, 


Frequency(GHz) 


FIG.  4.  (Color  online)  n(oj)  determined  from  the  focusing  experiments  using 
PhC  lens.  Note  that  the  theoretical  prediction  (solid  line)  matches  excel¬ 
lently  with  the  experimental  results. 


and  these  windows  move  downward  as  the  radius  of  curva¬ 
ture  is  increased.  The  inhomogeneity  of  the  PhC  lens  leads  to 
corrugation  on  the  concave  surface.  Windowing  of  refractive 
indices  for  each  radius  of  curvature  could  be  due  to  corruga¬ 
tion,  as  the  conservation  of  wave  vector  on  the  refracting 
surface  is  considerably  effected  from  such  corrugation  lead¬ 
ing  to  no  focus  point.  However,  in  the  regions  where  focus¬ 
ing  is  achieved,  corrugation  has  minimal  effect  on  the  sharp¬ 
ness  of  the  image. 

The  present  lens  has  several  advantages  when  compared 
to  the  one  with  positive  index.  Lenses  with  reduced  geomet¬ 
ric  aberrations  produce  sharper  image  with  enhanced  resolu¬ 
tion  and  find  numerous  applications.  For  any  value  of  n  <  0, 
the  radius  of  curvature  of  the  left-handed  PhC  lens  is  always 
larger  than  that  of  its  counterpart,  a  positive  index  planocon¬ 
vex  lens.  Larger  radius  of  curvature  gives  the  advantage  of 
reduced  aberration  in  the  image  formed.  Secondly,  a  PhC 
lens  having  the  same  focal  length  as  that  of  a  conventional 
lens  weighs  far  less,  and  is  suitable  for  space  applications. 
The  tailor  made  refractive  index  achievable  in  PhC 
materials12  allows  further  control  on  the  focal  length  and 
thereby  helps  reduce  the  length  of  the  optical  systems. 

Bandwidth  for  obtaining  a  sharp  focus  point  is  a  crucial 
parameter  that  decides  the  eventual  applications  of  the  left- 
handed  lenses.  The  present  PhC  lens  reveals  a  wide  band¬ 
width  of  2  GHz,  which  is  22.7%  at  the  current  operating 
frequencies.  In  comparison  to  a  planoconcave  lens  made  of 


wave.  Scale:  On  the  left  side  intensity  varies  from  0.005  to  0.055,  on  the 

right  side  the  real  part  of  S21  from  -0.025  to  0.015.  Dimensions  of  the  “  TlG.  5.  (Color  online)  Band  structure  of  the  PhC  calculated  from  plane 
figure  are  49  X  34  cm2.  wave  expansion  method. 
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the  MM,  the  PhC  bandwidth  is  much  larger.  Due  to  the  reso¬ 
nant  nature  of  the  MM  the  bandwidth  is  usually  restricted  to 
a  narrow  region  and  the  dispersion  is  stronger/3  The  weaker 
dispersion  in  the  PhC  makes  it  a  better  candidate  for  focus¬ 
ing  a  pulse  or  broadband  radiation. 

In  conclusion  the  feasibility  of  designing  a  broadband 
left-handed  lens  is  experimentally  demonstrated.  Focusing  of 
plane  waves  by  the  planoconcave  PhC  lens  is  achieved  for 
three  different  radii  of  curvature.  The  focal  length  follows 
the  standard  laws  of  geometrical  optics  applied  with  negative 
refraction.  Further,  the  measured  values  of  refractive  indices 
of  the  lens  are  in  complete  agreement  with  those  determined 
from  band  structure  calculations. 
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Remarkably  slow  propagation  of  microwaves  in  two  different  classes  of  left-handed  materials  (LHM’s)  is 
reported  from  microwave-pulse  and  continuous-wave  transmission  measurements.  Microwave  dispersion  in  a 
composite  LHM  made  of  split-ring  resonators  and  wire  strips  reveals  group  velocity  vg  ~  c/50,  where  c  is  the 
free-space  light  velocity.  Photonic  crystals  (PhC’s)  made  of  dielectric  A1203  rods  reveal  u4,~c/10.  Group 
delay  dispersion  of  both  the  composite  LHM  and  PhC’s  determined  from  the  experiment  is  in  complete 
agreement  with  that  obtained  from  theory.  The  slow  group  velocities  are  quantitatively  described  by  the  strong 
dispersion  observed  in  these  materials. 

DOI;  10. 1 103/PhysRevB.72.033 1 10  PACS  number(s):  41.20.Jb,  42.70.Qs,  72.30.+q 


Recent  observation  of  left-handed  electromagnetism  in  a 
composite  metamaterial1  (CMM)  made  of  split-ring  resona¬ 
tors  and  wire  strips  and  photonic  crystals  (PhC’s)  comprised 
of  periodic  arrangements  of  metallic  or  dielectric  elements 
has  revealed  the  unique  properties  of  these  left-handed  ma¬ 
terials  (LHM’s).  Negative  refraction1-4  and  imaging  by  flat 
lens,5'6  which  are  consequences  of  left-handed  electromagne¬ 
tism,  are  well  established  in  both  CMM  and  PhC’s.  The 
CMM  simultaneously  possesses  negative  permittivity  e'<0 
and  permeability  fi'  <  0  in  a  small  frequency  window,  result¬ 
ing  in  a  refractive  index  n'  < 0,  since  n~  \yl X  < 0  with 
n(a>)=n’  +  in",  +  and  e(w)=e'  +  re".  The  PhC, 

although  locally  has  positive  /c'(to)  and  e'(w),  possesses  ef¬ 
fective  negative  refractive  index  n'{u>)< 0  in  certain  fre¬ 
quency  ranges  due  to  its  dispersion  characteristics. 

The  electromagnetic  (EM)  waves  are  expected  to  experi¬ 
ence  delay  in  the  LHM  due  to  the  strong  dispersion,  first 
noted  by  Veselago.7  Another  peculiar  EM  property  of  the 
LHM  is  backward  wave  propagation  with  opposing  phase 
and  group  velocity  directions.  Fundamentally,  investigation 
of  group  velocity  is  interesting  as  the  wave  propagation  di¬ 
rection  and  speed  are  determined  by  the  group  velocity  rather 
than  the  phase  velocity.  Also,  low  group  velocities  would 
allow  control  of  light  propagation  and  could  lead  to  several 
applications  including  delay  line  filters  and  phase  shifters, 
from  microwave  to  optical  frequencies.  Photonic  dispersion 
determined  previously  for  multilayer  of  Fibanocci  quasi 
crystals8,9  revealed  a  slowing  down  of  the  wave  packets 
close  to  the  photonic  band  edge  with  vg  ~  c/5.  Slow  v„  of  the 
order  of  order  of  c/3  was  observed  in  Si/Si02  photonic 
crystals,10  and  similar  values  are  obtained  on  AlGaAs-based 
photonic  crystal  slabs.11  Besides,  a  recent  theoretical 
investigation12  has  proposed  subtuminal  propagation  in  the 
LHM. 

In  this  report  we  present  the  results  of  pulse  propagation 
in  two  different  classes  of  left-handed  materials:  the  CMM 
and  PhC’s.  The  results  show  that  the  group  velocity  is  as  low 
as  c/50  in  the  CMM  and  c/10  in  PhC’s.  The  experimental 
results  of  group  velocity  dispersion  are  in  good  agreement 
with  the  theory. 


The  CMM  was  made  by  interleaving  parallel  arrays  of 
split-ring  resonators  (SRR’s)  and  wire  strips  (WS’s)  etched 
out  on  a  circuit  board  using  photolithography.  The  circuit 
board  material  is  made  of  Taconic  FR-35  with  s' =4.7  in  the 
microwave  region.  The  cell  size  is  0.5  cm,  the  distance  be¬ 
tween  the  interleaved  arrays.  An  array  of  SRR’s  and  WS’s 
should  result  in  an  effective-negative-permeability  and 
negative-permittivity  material.  Therefore  a  lattice  of  such 
material  produces  a  negative-index  metamaterial  (NIM)  that 
has  effective  n'  < 0.  Further  details  of  the  fabrication  of  the 
CMM  are  given  in  Refs  1,  13,  and  15. 

The  microwave  PhC’s  consist  of  an  array  of  cylindrical 
A1203  rods,  having  dielectric  constant  e'=9,  radius 
0.316  cm,  and  height  1.25  cm  arranged  on  a  square  lattice. 
Two  different  PhC’s  having  different  r/a  ratios,  r  being  the 
radius  and  a  lattice  constant,  were  fabricated.  For  PhCl, 
r/a=0.175,  and  for  PhC2,  it  was  0.35.  These  left-handed 
PhC’s  manifest  tailor-made  refractive  index14  with  high 
dn/dco  which  is  well  established  from  continuous-wave  (cw) 
microwave  experiments. 15 

Figure  1  shows  a  schematic  diagram  of  the  experimental 
setup  used  for  the  pulse  measurements  carried  out  in  the  time 
domain.  The  setup  consists  of  an  HP  70820A  microwave 
transition  analyzer  (MTA)  connected  to  an  HP  8341 B  syn¬ 
thesized  sweeper  (SS).  A  Princeton  pulse  generator  FG100  is 
employed  to  get  a  pulse  of  width  35  ns.  The  pulse  thus  gen¬ 
erated  with  a  earner  frequency  in  the  microwave  region  is 
made  to  pass  through  two  transmission  lines  (A  and  B)  using 
a  splitter.  Two  isolators  are  connected  to  the  beam  splitter  to 
prevent  reflected  signal  interfering  with  the  signals  fed  to  the 
reference  and  sample  transmission  lines.  The  signals  are  ana¬ 
lyzed  by  the  two  MTA  input  channels.  To  enhance  the  sen¬ 
sitivity  of  the  measurement,  the  difference  in  the  time  delay 
between  the  reference  and  transmission  lines,  with  no  sample 
loaded,  is  set  to  zero.  An  X-band  rectangular  waveguide 
transmission  line  is  used  for  measurements  on  the  CMM. 
The  PhC  was  assembled  in  a  parallel  plate  waveguide  in 
which  the  propagation  mode  is  TM  with  electric  field  parallel 


to  the  AI1O3  rods.  In  all  the  measurements,  the  background 
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FIG.  1 .  Photographs  of  composite  metamaterial  (a)  and  photonic 
crystal  (b).  Schematic  diagram  of  the  pulse  delay  measurement 
setup  (c). 


time  delay  caused  by  the  transmission  lines  and  rectangular 
waveguide  is  well  accounted  for  to  obtain  the  time  delay  due 
to  the  sample  alone.  The  experiment  is  validated  by  measur¬ 
ing  pulse  delay  in  an  empty  waveguide,  which  matches  with 
the  predicted  value  as  shown  in  Fig.  2. 

Previous  cw  measurements'5  on  the  CMM  indicate  that 
the  LH  transmission  band  lies  between  9.8  and  10.4  GHz. 
Therefore  pulse  measurements  were  carried  out  in  this  region 
on  CMM  samples  of  two  different  lengths:  (=1.5  cm  (four 
slabs)  and  Z=1  cm  (three  slabs).  Figure  2  shows  the  group 
velocity  dispersion  in  the  CMM  of  two  different  lengths  de¬ 
termined  from  pulse  delay  using  vg-llt,  where  t  is  the  time 
delay.  A  similar  group  velocity  dispersion  in  both  samples 
indicates  that  there  is  no  dependence  on  the  sample  thick¬ 
ness.  The  vg  is  very  low  and  varies  from  0.08  to  0.01  be¬ 
tween  9.8  and  10.3  GHz.  vg  appears  to  diverge  near  the  low¬ 
est  frequencies.  These  frequencies  correspond  to  the  edge  of 
the  CMM  passband  in  which  the  left-handed  nature  is  not 
well  established.  Consequently  only  the  results  in  a  fre¬ 
quency  window  9.9- 10.3  GHz  are  considered  representative 
of  the  CMM. 

Pulse  delay  in  homogenous  and  inhomogeneous  material 
media  can  occur  due  to  scattering  of  waves.  Also,  at  the 
interface  of  vacuum  and  the  material  medium  a  pulse  suffers 
from  a  short  delay  due  to  the  transient  time  needed  for  the 
transmitted  wave  to  reorganize  and  eventually  pass  through 
the  PhC.  However,  the  time  delay  caused  by  these  two 
mechanisms  is  short  and  cannot  account  for  the  long  delay 
times  observed.  The  present  results  can  be  understood  from 
the  dispersion  characteristics  of  the  waves  in  the  CMM. 

Analytic  calculations  suggest  that  effective  s(«)  of  WS 
(Ref.  16)  and  of  SRR  (Ref.  17)  can  be  described  by 


e(w),/Z(w) : 


0)2  -  °>ej«a  +  io)7e,m 


where  co=27 r/,  a>e0  (wm0)  is  the  low-frequency  edge  of  the 
electric  (magnetic)  forbidden  bands,  v>ep  («m;))  is  electric 
(magnetic)  plasma  frequencies,  and  ye  (ym)  is  the  corre¬ 
sponding  damping  factor.  The  index  of  refraction  n'(co) <0 


FIG.  2.  Group  velocity  dispersion  of  CMM  in  a  waveguide. 
Slowest  Vg  of  c/50  is  obtained  at  10.3  GHz  in  this  configuration. 
Solid  line  is  a  fit  to  the  analytical  model. 


in  the  range  max((om0,a)eo)<w<min((omp,<ot,p).  Due  to  the 
resonance  effect  just  below  the  magnetic  and  electric  plasma 
frequencies,  both  permittivity  and  permeability  functions  un¬ 
dergo  large  changes  with  frequency,  which  results  in  unusual 
negative  values  for  these  two  material  parameters.  Conse¬ 
quently  a  negative  value  can  be  assigned  to  the  refractive 
index  «'(&>)  when  both  /x'(a>)  and  s'(o))  are  negative,  as  well 
as  a  strong  dispersion  results.  A  notable  feature  of  the  group 
velocity  dispersion  is  the  decrease  in  vg  with  the  increase  in 
frequency,  which  validates  the  nature  of  the  analytical  ex¬ 
pression  in  accordance  with  the  causality  requirement.  The 
present  vg  dispersion  is  robust  and  compares  well  with  that 
determined  from  the  cw  technique.15  The  excellent  agree¬ 
ment  between  these  two  values  is  evident  from  Fig.  2.  Also, 
a  fit  to  the  experimental  data  using  the  analytical  model  and 
vg = fil\k{n  +  codn/dw)]  is  shown  in  the  figure,  where  k 
=na)/c,  /3='Jk2-(v/a)1,  and  a =23  cm,  the  width  of  the 
x-band  waveguide.  It  can  seen  that  the  vg  determined  from 
the  model  follows  the  experimental  data  very  well  for  the  fit 
parameters  feo= 5  GHz,  fmo= 8  GHz,  fep= 25  GHz,  fmp 
=  10.9  GHz,  ye=L3  GHz,  and  y;„= 0.043  GHz.  These  values 
are  consistent  with  that  obtained  from  the  cw  measurement.15 
The  free-space  group  velocity  determined  using  the  model 
varies  from  0.052  to  0.034  (in  the  units  of  1/c)  in  the  region 
9.8-10.3  GHz. 

In  Fig.  3,  og  of  PhCl  of  thickness  10  layers  for  the  pulse 
propagation  along  the  TY  direction  of  the  first  Brillouin  zone 
is  shown.  The  slow  wave  propagation  with  vg  varying  form 
0.2  to  0.35  in  the  region  7.73-9.4  GHz  is  evident.  Figure  4 
shows  the  group  delay  dispersion  for  three  different  samples 
of  thickness  8,  12,  and  16  layers  of  PhC2.  For  all  the  three 
measurements  the  pulse  propagation  is  along  the  TM  direc¬ 
tion  of  the  first  Brillouin  zone.  Note  that  the  vg  is  very  low 
and  varies  from  0.45  to  0.1  between  7  and  8.3  GHz.  The 
results  also  indicate  that  there  is  no  dependence  of  vg  on 
sample  thickness.  As  in  the  case  of  PhCl,  PhC2  also  indi¬ 
cates  slow  group  velocities. 

The  low  vg  obtained  on  all  the  samples  of  PhCl,  PhC2, 
and  CMM  is  striking  and  forms  a  salient  feature  of  this  com¬ 
munication.  In  the  case  of  PhC’s,  the  slow  vg  can  be  attrib¬ 
uted  to  the  dispersion  characteristics  of  the  waves,  which  can 
be  determined  from  the  band  structure. 
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FIG.  3.  Group  velocity  dispersion  of  lefthanded  PhCl  with 
r/o=0. 175.  Note  the  low  vg  with  the  lowest  value  of  c/5  at 
9.3  GHz. 


FIG.  4.  Group  velocity  dispersion  of  right-handed  PhC2  with 
rla= 0.35.  Note  the  very  low  vg  with  the  lowest  value  of  c/10  at 
8.3  GHz. 


The  band  structures  of  the  PhCl  and  PhC2  are  shown  in 
Fig.  5.  We  have  investigated  the  band  structure  of  both  the 
PhC’s  for  left-handed  or  right-handed  behavior  of  micro- 
waves  following  Refs.  2,  4,  and  18.  The  strong  modulation 
in  the  PhC  would  result  in  guided  modes  with  waves  having 
different  phase  velocities  at  different  frequencies.  The 
phase  and  group  velocities  in  a  medium  are  given  by  vp 
=  (c/\nl,\)ij  with  kj=kp/\kp\  and  vg=Vg(o.  An  effective  refrac¬ 
tive  index  can  be  defined,  np=sgn(vg  •  kf)c\kj\  /  <w,  and  calcu¬ 
lated  from  the  band  structure.  The  sign  of  np  is  determined 
from  the  behavior  of  the  EFS, 2,19,20  which  are  drawn  from 
the  dispersion  curve  of  the  band  structure. 

In  the  case  of  PhCl  with  r/a=0.175  careful  observation 
of  the  band  structure  (solid  line  in  Fig.  5)  reveals  that  in  the 
second  band  near  the  band  edges  where  the  light  line  crosses 
the  bands,  the  group  velocity  is  antiparallel  to  the  propaga¬ 
tion  vector  k  and  the  refractive  index  is  negative.  This  is  a 
characteristic  of  the  left-handed  behavior  in  the  PhC.  The 
group  velocity  dispersion  determined  from  the  band  structure 
(solid  line)  is  shown  in  Fig.  3.  The  experimental  result  strik¬ 
ingly  matches  with  the  theory.  A  key  feature  of  the  band 
structure  is  the  flattening  of  the  bands  at  the  band  edges.  The 
vg  in  these  regions  is  expected  to  be  low.  Below  the 
band  edge  7.7  GHz,  vg=0  and  no  propagation  takes  place. 
It  is  interesting  to  note  that  as  frequency  decreases  from 
8.5  to  7.7  GHz,  vg  decreases  and  both  vg  and  the  refractive 
index  approach  zero.2  The  discrepancy  between  the  experi¬ 
mental  and  theoretical  vg  near  the  band  edge  7.7  GHz  is 
due  to  a  drastic  reduction  in  signal  strength  at  the  band 
edges,  which  results  in  larger  error  in  the  measurement 
of  vg. 

Following  a  similar  band  structure  analysis  for  PhC2  it 
can  be  deduced  that  in  the  region  7-8.5  GHz  the  electro¬ 
magnetism  is  righthanded.  In  this  region  the  directions  of  k 
and  vg  such  that  vg-kf>0.  Figure  4  shows  the  group  velocity 
dispersion  determined  from  band  structure.  Evidently,  as  in 
the  case  of  PhCl,  excellent  agreement  between  the  theoreti¬ 
cal  and  experimental  results  is  achieved  in  the  entire  fre¬ 
quency  range.  It  is  important  to  note  that  maximum  slowing 
down  of  the  vg  is  obtained  close  to  the  band  edge  8.3  GHz.  A 


similar  low  vg  value  for  PhCl  could  have  been  obtained  but 
is  prevented  by  the  very  weak  signal  strength. 

In  conclusion,  low  group  velocities  are  observed  in  left- 
handed  metamaterial  and  photonic  crystals.  These  results 
validate  the  behavior  of  the  waves  expected  from  the  disper¬ 
sion  characteristics  in  these  media.  A  quantitative  description 
of  the  experimental  results  is  obtained  using  the  analytical 
model  in  the  case  of  the  CMM  and  band  structure  in  the  case 
of  PhC’s.  An  important  feature  is  a  strong  slowing  down  of 
the  waves  at  the  band  edges  of  the  PhC’s.  The  group  velocity 
reduction  observed  in  the  photonic  crystals  is  more  than  3 
times  larger  than  that  of  Fibanocci  quasicrystals  and  9  times 
that  of  colloidal  photonic  crystals  of  polystyrene.21  Such  low 
vg  are  due  to  strong  modulation  (with  a  dielectric  contrast 
9-1)  in  the  present  PhC’s  by  way  of  weak  modulation  in  the 
colloidal  polystyrene  (contrast  2.5-1)  PhC.  The  present  in¬ 
vestigation  indicates  that  very  slow  group  velocities,  lower 
than  that  presented  here,  are  possible  in  metamaterials  and 
photonic  crystals  having  low  transmission  loss.  These  mate¬ 
rials  can  be  tailor  made  to  have  a  strong  dispersion  with  very 
high  dn/dto  to  result  in  lower  group  velocities.  The  availabil¬ 
ity  of  high-permittivity  materials  with  minimal  losses  facili¬ 
tates  the  design  of  tailor  made  photonic  crystals  with  slow 
vg.  The  ability  to  mold  and  slow  light  leads  to  several  appli- 


FIG.  5.  TM-mode  band  structure  of  square  lattices  PhCl  and 
PhC2.  Note  the  flattening  of  the  bands  at  the  band  edges  which 
results  in  slow  group  velocities. 
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cations  of  metamaterials  and  photonic  crystals  which  are  not 
possible  with  conventional  materials.  Important  of  them  are 
delay  line  filters  and  phase  shifters.  The  slow  group  veloci¬ 
ties  obtained  at  microwave  frequencies  can  also  be  obtained 
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at  optical  frequencies  by  scaling  the  material  parameters  to 
micron  sizes. 
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We  demonstrate  the  negative  refraction  of  microwaves  in  a  metallic  photonic  crystal  prism.  The 
spectral  response  of  the  photonic  crystal  prism,  which  manifests  both  positive  and  negative  refraction, 
is  in  complete  agreement  with  band-structure  calculations  and  numerical  simulations.  The  validity  of 
Snell’s  law  with  a  negative  refractive  index  is  confirmed  experimentally  and  theoretically.  The  negative 
refraction  observed  corresponds  to  left-handed  electromagnetism  that  arises  due  to  the  dispersion 
characteristics  of  waves  in  a  periodic  medium.  This  mechanism  for  negative  refraction  is  different  from 
that  in  metamaterials. 
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The  optical  properties  of  materials  that  are  transparent 
to  electromagnetic  (EM)  waves  can  be  characterized  by 
an  index  of  refraction  n.  Given  the  direction  of  the 
incident  beam  Bx  at  the  interface  of  vacuum  and  the 
material,  the  direction  02  of  the  outgoing  beam  can  be 
determined  using  Snell’s  law  sin#|  =  n  sin02.  All  natu¬ 
rally  available  materials  possess  a  positive  refractive 
index  n  >  0.  It  was  observed  recently  that  in  certain 
composite  metamaterials  EM  waves  bend  negatively  [1] 
and,  consequently,  a  negative  index  of  refraction  n  <  0 
can  be  assigned  to  such  materials  without  violating 
Maxwell’s  equations.  For  homogeneous  media,  a  negative 
refraction  (02  <0,  0X>  0)  necessarily  requires  that 
n  <  0.  This  negative  bending  allows  considerable  control 
over  light  propagation  and  opens  the  door  for  new  ap¬ 
proaches  to  a  variety  of  applications  from  microwave  to 
optical  frequencies. 

Negative  index  media  exhibit  some  unusual  propaga¬ 
tion  characteristics  of  the  EM  waves.  The  most  striking 
property  is  that  of  left-handed  electromagnetism  (LHE) 
since  the  electromagnetic  fields  E  and  H  and  the  wave 
vector  k  form  a  left-handed  triplet  Consequently,  the 
energy  flow  represented  by  the  Poynting  vector  S  =  EX 
H  is  antiparallel  to  the  wave  vector  k,  so  that  S  •  k  <  0.  In 
contrast,  for  conventional  n  >  0  materials,  E,  H,  k  form  a 
right-handed  triplet  corresponding  to  right-handed  elec¬ 
tromagnetism  (RHE),  and  S  ■  k  >  0. 

A  material  possessing  simultaneously  negative  permit¬ 
tivity  e  <  0  and  permeability  jjt<  0  can  be  shown  to 
necessarily  have  n  =  <  0  [2,3].  Recently,  negative 

refraction  was  demonstrated  in  a  quasihomogeneous 
metamaterial  [1,4,5]  with  unit  cell  dimensions  less  than 
the  wavelength,  consisting  of  interleaved  arrays  of  wires 
(e  <  0)  and  split  ring  resonators  ( fi  <  0).  However,  these 
materials  are  highly  absorptive,  and  unlikely  to  be  scaled 
to  three  dimensions  or  to  smaller  sizes  suitable  for  appli¬ 
cations  at  optical  frequencies  [1,6,7]. 


It  has  since  been  proposed  that  negative  refraction 
can  be  achieved  in  photonic  crystals  (PC),  which  are 
inhomogeneous  periodic  media  with  a  lattice  constant 
comparable  to  the  wavelength  [8,9],  A  PC  is  an  artificial 
structure,  usually  made  of  a  dielectric  or  metal,  designed 
to  control  photons  similar  to  the  way  a  solid  state 
crystal  controls  electrons.  Locally  both  e,  fi  >  0  every¬ 
where  in  the  PC.  The  physical  principles  that  allow 
negative  refraction  in  the  PC  arise  from  the  disper¬ 
sion  characteristics  of  wave  propagation  in  a  periodic 
medium,  and  are  very  different  from  that  of  the 
metamaterial  in  Ref.  [1],  In  this  Letter,  we  present 
experimental  evidence  of  negative  refraction  of  micro- 
waves  in  a  metallic  photonic  crystal.  Parallel  theoretical 
and  numerical  investigations  of  the  band  structure 
and  simulations  of  wave  propagation  through  the  PC 
prism  are  compared  with  the  experimental  results. 
Exceptionally  good  agreement  is  found  between  the  ex¬ 
periments,  band-structure  calculations,  and  wave  refrac¬ 
tion  simulations. 

The  microwave  photonic  crystal  fabricated  in  the 
present  work  is  an  array  of  cylindrical  copper  rods  of 
height  1.26  cm  and  radius  0.63  cm  forming  a  triangular 
lattice.  The  ratio  of  the  radius  r  to  lattice  constant  a  is 
r/a  =  0.2.  Refraction  measurements  are  carried  out  in  a 
parallel  plate  waveguide  made  of  a  pair  of  metallic  plates. 
The  excitation  in  the  parallel  plate  waveguide  is  a  trans¬ 
verse  magnetic  (TM)  mode  up  to  12  GHz  such  that  the 
electric  field  E  is  parallel  to  the  rod  axis.  Microwave 
absorbers  are  placed  on  both  sides  to  avoid  spurious 
reflections  and  to  collimate  the  propagating  beam  of 
width  9  cm,  which  is  incident  normally  on  a  right  angle 
prism  of  PC  [see  Fig.  1(a)].  On  the  far  side,  a  dipole 
antenna  attached  to  an  X-Y  robot  maps  the  electric  field 
E.  Accurate  angles  of  refraction  are  obtained  by  fitting 
the  emerging  wave  front  with  a  plane  wave,  with  the 
refraction  angle  d2  as  a  fit  parameter. 
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FIG.  1  (color),  (a)  Schematic  diagram  of  microwave  refraction  experimental  setup  (not  to  scale),  (b)  Propagation  vectors  for 
positive  and  negative  refraction,  (c)-(f)  Microwave  electric  field  maps  in  the  far  field  region,  (c)  Negative  and  (e)  positive  refraction 
by  the  metallic  PC  prism  for  the  incident  beam  along  T  — » JC  forming  an  angle  of  incidence  30°  with  the  refraction  surface.  WF  = 
wave  front  with  respect  to  refracting  surface,  (d)  Negative  refection  for  the  incident  beam  along  T  — ►  M  forming  an  angle  of 
incidence  60°.  (f)  Positive  refraction  by  a  polystyrene  prism.  In  all  the  field  maps,  approximate  area  of  each  field  map  is  43  X 
40  cm2;  the  PC  prisms  and  incident  beams  are  schematic  and  do  not  correspond  to  the  actual  sizes  used. 


The  refraction  experiment  is  validated  by  the  data  of  a 
polystyrene  prism  having  similar  dimensions  as  that  of 
the  PC  prisms.  In  Fig.  1(f),  the  direction  of  the  emerging 
beam  can  be  clearly  seen  at  an  angle  02  =  +52.2°  from 
the  normal  to  the  surface  of  refraction,  corresponding  to  a 
positive  refractive  index  n  =  1.58(nsin0i  =  sind2)  for  an 
incident  angle  0X  =  30°. 

The  measurements  on  the  triangular  lattice  PC  were 
carried  out  with  the  incident  wave  vector  k;  along  direc¬ 
tions  T  — ►  K  and  T  — »  M  of  the  first  Brillouin  zone  (these 
directions  in  real  space  and  in  reciprocal  wave-vector 
space  are  shown  in  Figs.  1(b)  and  3).  The  angles  of 
incidence,  0lK  =  30°  for  T  — ►  K  and  0XM  —  60°  for  T  — ► 
M,  are  chosen  in  such  a  way  that  the  periodicity  on  the 
surface  of  refraction  is  minimum  to  prevent  higher  order 
Bragg  diffraction  [5].  Figure  1(c)  illustrates  the  negatively 
refracted  wave  front  with  an  angle  of  refraction  02K  = 
11.5°  for  /  =  9.77  GHz  with  incidence  along  T  — ►  A". 
Using  Snell’s  law,  neff  sin#1A-  =  sin02/r  with  0XK  =  30° , 
02K  =  -11.5°,  we  obtain  an  effective  refractive  index  of 
neff  =  —0.4  at  this  frequency.  A  second  wave  front  can 
also  be  seen  emerging  from  the  top  edge  of  the  PC.  We 
attribute  this  wave  front  to  the  edge  effect  due  to  the  finite 
sample  size.  The  negative  refraction  reported  in  the 
present  photonic  crystal  has  been  demonstrated  for  the 
first  time  in  an  inhomogeneous  system  [10].  Figure  1(e) 
shows  positive  refraction  at  /  =  6.62  GHz.  In  Fig.  1(d), 
an  illustrative  example  of  negative  refraction  at  10.4  GHz 
for  F  — ►  M  (0lM  —  60°)  is  shown. 
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An  understanding  of  negative  refraction  and  its  relation 
to  the  left-handed  behavior  of  electromagnetic  waves  in  a 
photonic  crystal  can  be  achieved  by  examining  the  band 
structure  of  an  infinite  PC.  We  have  calculated  the  band 
structure  employing  standard  plane  wave  expansion 
methods  using  over  2000  plane  waves  [11].  The  2D  band 
structure  for  a  triangular  lattice  PC  with  r/ a  —  0.2  is 
shown  in  Fig.  2(a). 

For  a  plane  wave  with  wave  vector  k\  and  frequency  o> 
incident  normal  to  an  air-PC  interface,  the  wave  vector  kf 
inside  the  PC  is  parallel  or  antiparallel  to  kx  as  deter¬ 
mined  by  the  band  structure.  If  d<o/d\kf\ >  0,  kf  is 
parallel  to  kx  and,  consequently,  the  EM  field  in  the  PC 
is  right  handed  (RHE).  Otherwise,  kf  is  antiparallel  to  kx 
and  the  EM  field  in  the  PC  is  left  handed  (LHE).  For  a 
general  case,  the  phase  and  group  velocities  in  a  medium 
are  vp  =  (c/\hp\)kf  with  kf  —  kf/\kf\  and  vg  =  V^tw.  It 
can  be  proved  analytically  that  the  direction  of  group 
velocity  vg  in  an  infinite  PC  coincides  with  that  of  the 
energy  flow  S  [12],  An  effective  refractive  index  can  be 
defined  np  =  sgn(ufi  •  £/)[(c|£/|)/<w]  and  calculated  from 
the  band  structure  [13].  The  sign  of  np  is  determined  from 
the  behavior  of  the  equifrequency  surfaces  (EFS).  EFS 
plots  for  the  first  and  second  bands  of  the  triangular 
lattice  are  shown  in  Fig.  3(a).  The  center  in  the  plots 
corresponds  to  the  center  of  the  first  Brillouin  zone. 
The  EFS  that  move  outwards  from  the  center  with 
increasing  frequency  correspond  to  the  RHE  with  vg  • 
kf  >  0  and  inward  moving  surfaces  correspond  to  the 
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FIG.  2  (color  online),  (a)  Band  structure  computed  for  the 
triangular  lattice  metallic  PC  with  r/a  =  0.2.  (b)  Microwave 
transmission  amplitude  |S2il  vs  frequency  /  (GHz)  on  the  far 
side.  Note  the  band  cutoff  below  6.2  GHz,  transmission  be¬ 
tween  5.8-11.1  GHz,  and  band  gap  region  11.1-11.3  GHz  are  in 
excellent  agreement  with  Fig.  2(a).  (c)  Refractive  index  np 
determined  from  the  experimental  results  for  a  beam  incident 
along  T  — ►  K  (A)  and  T  — *  M  (*).  For  np<  0  (np  >  0),  the 
electromagnetism  is  left  handed  (right  handed),  and  the  refrac¬ 
tion  angle  is  negative  (positive).  The  close  match  between  the 
experiment  (A  and  *)  and  theory  (dashed  and  solid  lines)  is 
evident  and  shows  agreement  with  Snell’s  law. 


LHE  with  vs,-kf<  0.  In  the  case  of  LHE  (RHE)  con¬ 
servation  of  kf  component  along  the  surface  of  refraction 
would  result  in  negative  (positive)  refraction.  The  result¬ 
ing  refractive  index  np  determined  from  the  band  struc¬ 
ture  and  EFS  for  a  beam  incident  along  both  T  — ►  K 
(dashed  line)  and  Y  —*  M  (solid  line)  is  shown  in 
Fig.  2(c).  Note  that  negative  refraction  is  predicted  for 
regions  in  the  second  and  third  bands  and  positive  refrac¬ 
tion  in  the  first  and  fourth  bands. 

In  the  following,  we  describe  the  salient  features  of  the 
experimental  results  and  comparison  to  band  structure. 

(i)  In  the  first  band  between  5.8-7.7  GHz,  the  EFS  move 
outward  with  increasing  frequency,  so  that  np>  0  corre¬ 
sponding  to  RHE  with  Vg'kf>  0  (i.e„  vg  parallel  to  kf). 
The  representative  field  map  at  /  =  6.62  GHz  in  Fig.  1(e) 
confirms  the  positive  refraction  expected.  The  measured 
nef f  are  in  good  agreement  with  the  theoretical  calcula¬ 
tions  [Fig.  2(c)]. 

(ii)  In  the  second  band  between  7.7-11  GHz,  the  EFS 
move  inward  with  increasing  frequency,  consistent  with 
np<  0  corresponding  to  LHE  with  vg  •  kf  <  0  (vg  anti¬ 
parallel  to  kf).  The  illustrative  field  map  in  Fig.  1(c)  at 
/  =  9.77  GHz  shows  the  emerging  wave  front  in  the 
negative  direction.  The  experimental  results  of  refraction 
for  the  incident  beam  along  both  T  — *  K  [A  in  Fig.  2(c)] 
and  r  — *  M  (*)  are  in  excellent  agreement  with  the  band- 
structure  calculations  [dashed  and  solid  lines  in  Fig.  2(c)]. 

(iii)  In  certain  frequency  ranges  in  which  the  EFS  is 
circular  and  frequency  is  not  so  high,  the  first  order  Bragg 
diffraction  is  very  weak.  For  the  zeroth  order  Bragg  peak, 
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FIG.  3  (color),  (a)  The  EFS  for  the  first  and  second  bands.  The 
centers  of  the  hexagons  represent  the  centers  of  the  first 
Brillouin  zones  in  the  respective  bands.  Blue  (red)  color  rep¬ 
resents  lower  (higher)  frequency.  (b),(c)  Simulations  of  wave 
refraction  showing  the  wave  front  emerging  from  a  metallic 
PC:  (b)  positive  refraction  at  6.6  GHz  and  (c)  negative  refrac¬ 
tion  at  9.7  GHz.  The  PC  used  in  the  simulations  has  the  same 
size  as  that  used  in  the  experiment.  The  electric  field  is  plotted 
as  ElJ 3  for  better  visibility.  The  agreement  with  Figs.  1(e)  and 
1(c)  is  evident 


np  is  angle  independent  and,  consequently,  Snell’s  law  is 
applicable.  The  index  of  refraction  np  [see  Fig.  2(c)]  in  the 
region  8.8-11  GHz  determined  from  the  experimental 
wave  field  scans  for  different  angles  of  incidence, 
viz.  0iK  =  30°  for  Y  — ►  K  and  01M  =  60°  for  T  — *  M,  is 
nearly  angle  independent  due  to  the  circular  nature  of  the 
EFS,  confirming  the  validity  of  Snell’s  law  in  this  fre¬ 
quency  region.  Thus,  we  have  confirmed  the  validity  of 
Snell’s  law  both  experimentally  and  theoretically 
[Fig.  2(c)],  A  noteworthy  point  is  that  the  strong  contrast 
in  the  metallic  PC  leads  to  near  circular  EFS  [Fig.  3(a)] 
and  thus  results  in  negative  refraction  in  a  wider  fre¬ 
quency  range  than  that  of  a  dielectric  PC  [6,8,9], 

(iv)  The  band  cutoff  at  5.8  GHz  ,  transmission  between 
5.8-11.1  GHz,  and  band  gap  region  between  11.1— 
11.3  GHz,  all  of  which  are  observed  in  the  transmitted 
spectrum  shown  in  Fig.  2(b),  are  in  excellent  agreement 
with  the  band-structure  calculations  shown  in  Fig.  2(a). 

Direct  numerical  simulations  (see  Fig.  3)  of  wave  re¬ 
fraction  were  also  carried  Out  and  are  in  good  agreement 
with  the  experimental  results  and  band-structure  calcu¬ 
lations.  The  simulation  is  done  using  a  Green’s  function 
boundary  wall  approach  originally  developed  for  hard 
wall  potentials  in  quantum  mechanics  [14].  As  shown 
in  Fig.  3,  the  simulation  results  in  negative  refraction 
for  /  =  9.7  GHz  in  the  second  band,  and  positive  refrac¬ 
tion  for  /  =  6.6  GHz  in  the  first  band,  both  in  agreement 
with  the  experiment  (Fig.  1). 

127401-3 


33 


Volume  92,  Number  12 


PHYSICAL  REVIEW  LETTERS 


week  ending 
26  MARCH  2004 


In  conclusion,  we  have  experimentally  demonstrated 
negative  refraction  in  a  new  class  of  material,  metallic 
photonic  crystal.  A  major  feature  of  the  present  work 
is  the  extraordinary  level  of  control  exemplified  by 
the  convergence  between  the  experimental  data,  band- 
structure  calculations,  and  simulations.  This  means  that 
a  variety  of  tailor-made  structures  are  feasible  that  can  be 
designed  and  constructed.  There  are  numerous  possibil¬ 
ities  opened  up  by  the  present  results.  For  many  applica¬ 
tions  such  as  imaging,  one  requires  index  matching 
between  the  negative  index  material  and  surroundings 
(relative  index  n  —  ~  1)  accompanied  by  negligible 
losses.  These  requirements  are  more  easily  met  with  PC 
than  with  composite  negative  index  metamaterials. 
The  metallic  PC  offers  the  additional  advantages  of 
highest  dielectric  contrast,  low  attenuation,  and  the 
possibility  of  focusing,  which  are  evident  from  the 
present  data  Furthermore  the  microwave  PC  can  be 
easily  scaled  to  three  dimensions  [15],  and  to  optical 
frequencies  [16],  which  is  highly  unlikely  with  com¬ 
posite  negative  index  metamaterials  [1,6].  Thus,  the 
advantages  of  negative  refraction  and  left-handed  elec¬ 
tromagnetism  that  have  been  proposed  recently,  such  as 
imaging  by  flat  lenses  [17],  beam  steerers,  couplers,  and 
others,  as  well  as  some  entirely  new  possibilities,  are 
feasible  with  photonic  crystals  from  microwave  to  optical 
frequencies. 
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Negative  refraction  and  left-handed  electromagnetism  in  a  metallic  photonic  crystal  (PC)  wedge  are 
demonstrated  in  free  space  for  both  transverse  magnetic  and  electric  mode  propagation.  The 
experimental  results  are  in  excellent  agreement  with  numerical  calculations  based  on  the  band 
structure  with  no  fit  parameters  used  in  modeling.  The  results  demonstrate  precision  control  on  the 
design  and  fabrication  of  the  PC  to  achieve  tailor-made  refractive  indices  between  -0.6  and  +1.  © 

2004  American  Institute  of  Physics.  [DOI:  10.1063/1.1787959] 


Left-handed  electromagnetism  and  negative  refraction 
have  recently  been  observed  in  the  microwave  frequency 
range  in  composite  metamaterials  1  made  of  split  ring  reso¬ 
nators  and  wire  strips,  and  in  photonic  crystals  (PCs).  2,3 
Negative  refraction  in  left-handed  metamaterials  (LHM), 
when  used  effectively,  opens  the  door  for  approaches  to  a 
variety  of  applications  from  microwave  to  optical  frequen¬ 
cies.  An  interesting  application  of  negative  refraction  is  the 
superlensing  4  effect  by  a  flat  lens  with  no  curved  surfaces 
that  can  potentially  overcome  the  diffraction  limit  imposed 
by  conventional  lenses.  Indeed  a  flat  lens  without  optical  axis 
has  been  fabricated  5  recently  using  a  photonic  crystal  struc¬ 
ture.  However,  for  device  applications  it  is  important  to  have 
control  over  material  parameters  to  be  able  to  design  and 
predict  material  properties. 

Both  negative  refraction  and  left-handed  electromagne¬ 
tism  have  been  demonstrated  in  photonic  crystals  using  a 
parallel  plate  waveguide.6,7  The  two-dimensional  (2D)  paral¬ 
lel  plate  structures  used  in  all  previous  experiments  to  dem¬ 
onstrate  negative  refraction  confine  the  left-handed  material 
and  lead  to  spurious  edge  effects,  as  observed  in  Ref.  6. 
Close  to  the  surface  of  refraction,  the  emerging  wave  fronts 
interfere,  and  a  clear  negatively  refracted  beam  is  difficult  to 
observe.  Also,  the  coupling  of  the  photonic  crystal  modes 
with  the  incident  wave  is  different  for  transverse  magnetic 
(TM)  and  electric  (TE)  mode  propagations  and  is  not  well 
achieved  in  parallel  plates.  In  addition,  bandwidth  is  a  crucial 
element  for  device  applications  in  a  wide  frequency  range.  It 
is  therefore  essential  to  investigate  left-handed  electromag¬ 
netism  and  negative  refraction  in  a  photonic  crystal  sus¬ 
pended  in  free  space  for  both  TM  and  TE  mode  propaga¬ 
tions. 

In  this  letter  we  report  negative  refraction  for  both  TM 
(Ell  to  the  rod  axis)  and  TE  (El  to  the  rod  axis)  mode  propa¬ 
gation,  in  a  metallic  PC  prism  suspended  in  free  space.  Re¬ 
sults  show  that  a  PC  can  exhibit  negative  refraction  with 
tailor-made  refractive  indices  in  a  large  frequency  range.  The 
propagation  in  different  bands  of  the  PC  can  be  tuned  with 
frequency  to  obtain  either  negative  or  positive  refraction. 
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Thus,  the  present  tailor-made  PC  can  be  utilized  for  a  variety 
of  applications. 

The  microwave  photonic  crystal  consists  of  an  array  of 
cylindrical  copper  tubes  of  height  60  cm  and  outer  radius 
0.63  cm  arranged  on  a  triangular  lattice.  The  ratio  of  the 
radius  r  to  the  lattice  constant  a  was  r/a= 0.2.  Refraction 
experiments  were  performed  in  an  anechoic  chamber  of  di¬ 
mensions  5X8X4  m3  to  prevent  reflections  from  the 
walls.  A  square  X-band  horn  placed  at  3  m  from  the  PC 
acts  as  a  plane-wave  source  (Fig.  1).  Placing  a  piece  of  mi¬ 
crowave  absorber  with  a  6X6  in2,  aperture  in  front  of  the 
PC  narrows  the  incident  beam.  On  the  far  side  another  square 
horn  attached  to  a  goniometer  swings  around  in  two-degree 
steps  to  receive  the  emerging  beam.  Refraction  is  considered 
positive  (negative)  if  the  emerging  signal  is  received  to  the 
right  (left)  of  the  normal  to  the  surface  of  refraction  of  the 
PC.  Measurements  were  carried  out  with  the  incident  wave 
vector  kj  along  T— ►  Af(0, 1)  direction  of  the  first  Brillouin 
zone  of  the  PC  and  in  both  TM  and  TE  modes.  The  angle  of 
incidence  0-60°  for  T— is  chosen  in  order  to  minimize 
surface  periodicity  along  the  surface  of  refraction,  thus  elimi¬ 
nating  higher  order  Bragg  waves. 

Figure  2(a)  shows  a  plot  of  the  transmitted  intensity 
measured  at  different  angles  and  incident  frequencies  for  the 
TM  mode  propagation.  As  can  be  seen  lfom  the  figure,  be¬ 
tween  6  and  7.1  GHz  the  signal  emerges  on  the  positive  side 
of  the  normal  to  the  surface  corresponding  to  positive  refrac¬ 
tion.  No  transmission  is  observed  between  7.1  and  8.3  GHz. 
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FIG.  1 .  (Color  online).  Microwave  free-space  refraction  experiment  set  up 
in  an  anechoic  chamber.  Negative  or  positive  refraction  is  determined  from 
the  direcdon  of  the  emerging  signal  with  the  normal  to  the  surface  of 
refraction. 
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Refractive  index 

FIG.  2.  (Color  online),  (a)  Plot  of  refracted  wave  intensity  measured  at 
various  angles,  for  TM  mode  propagation,  (b)  Similar  plot  for  TE  mode. 
Note  the  negative  and  positive  refraction  observed  in  different  bands. 
Circles  are  the  theoretically  calculated  refractive  indices  corresponding  to  a 
zero  order  Bragg  wave,  and  stars  are  first-order  Bragg  wave,  both  of  which 
match  strikingly  with  the  experimental  results  without  any  fitting 
parameters. 

Above  8.3  upto  11  GHz,  two  signals  are  observed  on  the 
positive  and  negative  sides  of  the  normal.  The  negatively 
refracted  signal  is  strongest  around  10.7  GHz  and  the  posi¬ 
tively  refracted  signal  around  8.6  GHz.  Although  both  posi¬ 
tively  and  negatively  refracted  signals  are  observed,  with  the 
increase  in  frequency,  positive  signal  gets  weaker  while 
negative  signal  gets  stronger.  The  experimental  refractive  in¬ 
dex  n  is  obtained  from  applying  Snell’s  law  n 
=sin(0r)/sin(0,)  to  each  emerging  beam.  The  validity  of 
Snell’s  law  has  been  established  earlier  in  metallic  PCs.6 

We  have  also  carried  out  measurements  of  refraction  for 
TE  mode  propagation.  The  results  for  this  mode  are  shown 
in  Fig.  2(b).  Here,  negative  refraction  is  observed  between 
6.4-9. 8  GHz  and  positive  refraction  between  9.8-12  GHz. 

It  is  important  to  note  that  negative  refraction  is  possible  for 
both  TM  and  TE  modes;  such  a  freedom  in  the  choice  of 
modes  provides  a  crucial  advantage  of  using  the  metallic  PC 
over  the  split  ring  wire  array  metamaterial. 

The  physical  principle  behind  the  present  results  can  be 
understood  from  the  band  structure  of  the  metallic  PC.  We 
have  calculated  the  band  structure  of  the  triangular  lattice  PC 
employing  a  standard  plane-wave  expansion  method  using 
2500  plane  waves.  The  2D  band  structures  of  both  TM  and 

TE  are  shown  in  Fig.  3.  For  a  plane  wave  with  incident  wave 
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FIG.  3.  (Color  online).  Band  structures  of  the  triangular  metallic  PC  for  both 
TM  (dashed  lines)  and  TE  (solid  lines)  propagation  modes  are  presented. 

vector  k[  and  frequency  to  incident  normally  on  an  air-PC 
interface,  the  wave  vector  kf  inside  the  PC  is  parallel  or 
anti-parallel  to  kt  as  determined  by  the  band  structure  (Fig. 
4).  For  general  case,  the  phase  and  group  velocities  in  a 
medium  are  vp={d\np\)kf,  with  kj=kj/\kj\  and  vg=V^(». 
Note  that  the  direction  of  the  group  velocity  vg  in  an  infinite 
PC  coincides  with  that  of  the  energy  flow.  Conservation  of 
the  kf  component  along  the  surface  of  refraction  would  result 
in  positive  or  negative  refraction,  depending  on  whether  kf  is 
parallel  or  antiparallel  to  group  velocity. 

The  emerging  beam  can  be  written  as  'irl=a0e,kt>yr 
+ale'k'l'r,  where  kl0  and  ktX  represent  the  refracted  wave 
vectors  corresponding  to  the  zero-  and  first-order  Bragg 
wave  vectors  of  the  field  inside  the  PC.  Let  the  component 
of  k,  along  the  normal  to  the  surface  of  refraction  of  the 
prism  be  ktl  .  One  has  kl01  =  \laj2/c2-kj  sin2  8  and 

ktll  =  \jo>2/c2-(2Tr/a-kf  sin  0)2,  where  0=tt/3  and  a  is  the 
lattice  periodicity.  The  parallel  component  is  ka=kf  sin  6. 
The  refracted  angle  can  be  obtained  from  0r=tan_1(ifc,||//;,x) 
for  each  beam. 

From  the  band  structure  and  the  equifrequency  surfaces 
(EFS)  for  the  TM  mode  of  propagation,  negative  refraction  is 
predicted  for  the  second  and  third  band  regions,  with  positive 
refraction  in  the  first  band.  In  the  first  band  between 
6-7.1  GHz  the  EFS  move  outward  with  increasing 
frequency,  so  that  vg-kf> 0.  In  the  second  band  between 
8.3-11  GHz,  the  EFS  move  inward  with  increasing 
frequency,  consistent  with  np< 0,  corresponding  to  vg-kf<0 
(vg  antiparallel  to  kf).  The  band  gap  is  in  the  frequency  range 
7. 1-8.3  GHz  between  the  first  and  second  pass  bands  and 
from  11-11.2  GHz  between  the  second  and  third  bands. 


FIG.  4.  (Color  online.).  Directions  of  incident  transmitted  and  refracted 
3  (wave  vectors,  and  group  velocity  inside  the  PC.  Positive  or  negative  refrac¬ 
tion  can  be  observed  depending  on  whether  ug-iy> 0  or  v.-kj <0. 
to  AIP  license  or  copyright,  see  http://apl.aip.org/apl/copyright.jsp 
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An  effective  refractive  index  can  be  defined  np 
=&gn{vg-kf)c\kf\l <a)  and  calculated  from  the  band  structure. 
The  sign  of  np  is  determined  from  the  behavior  of  the  EFS. 
In  Figs.  2(a)  and  2(b)  the  refractive  indices  for  various  bands 
determined  from  the  theory  are  plotted.  The  close  match  be¬ 
tween  the  theoretical  predictions  and  experimental  data  is 
striking.  All  the  predicted  features  including  band  gap  and 
negative  and  positive  refraction  are  observed  in  the  experi¬ 
ment.  The  degeneracy  observed  in  the  case  of  TE  mode  is 
due  to  multiple  bands  for  a  single  frequency,  which  results  in 
multiple  propagation  k  vectors.  It  is  an  interesting  observa¬ 
tion  that  different  intensities  are  associated  with  different 
propagation  vectors  for  TE  mode.  For  this  mode  and  for 
Y->M  propagation  from  Fig.  3,  it  can  be  seen  that  the  higher 
the  slope  of  the  curve,  the  more  intense  the  beam. 

A  particular  feature  of  interest  is  the  bandwidth  for  nega¬ 
tive  refraction  and  lefthanded  electromagnetism.  From  Figs. 
2(a)  and  2(b)  it  can  be  deduced  that  the  bandwidth  for  the  TE 
mode  is  42%  and  for  TM  mode  27%.  In  comparison,  a  rela¬ 
tively  weakly  modulated  dielectric  PC  has  a  bandwidth  esti¬ 
mated  to  be  6%  (Ref.  3)  which  is  very  narrow  and  the  ex¬ 
perimentally  obtained  bandwidth  for  metamaterials  1,8  to  date 
is  only  10%  (Ref.  8).  The  present  bandwidths  for  both  TE 
and  TM  modes  are  higher  than  that  in  metamaterial.  Band¬ 
width  puts  stringent  restrictions  on  the  tenability  and  func¬ 
tional  range  of  the  devices  based  on  the  LHM.  In  particular, 
in  our  recent  work  we  have  shown  that  in  the  LHM,  electro¬ 
magnetic  (EM)  wave  propagation  is  slow  with  group  veloc¬ 
ity  of  0.02c,  9,10  where  c  is  the  velocity  of  the  EM  wave  in 
vacuum.  This  slow  group  velocity  combined  with  large  band¬ 
width  can  be  used  for  designing  a  delay  line  filter  with  a 
large  pass  band. 

In  the  case  of  the  TE  mode,  for  an  angle  of  incidence  of 
60°,  the  refractive  index  is  found  to  vary  from  0  to  -0.48, 
which  is  a  200%  change  for  a  frequency  change  of  42%. 


Such  a  large  dn/d(o  results  in  a  large  d<f>ldoi,  which  can  be 
used  in  designing  ultrasensitive  phase  shifters. 

In  conclusion,  negative  refraction  is  demonstrated  for 
both  TM  and  TE  mode  propagation  in  a  metallic  PC.  Almost 
400%  improvement  in  the  bandwidth  for  negative  refraction 
is  obtained  for  TE  mode  propagation  in  this  tailor-made  PC. 
The  strong  dispersion  observed  results  in  a  change  of  200% 
in  the  negative  refractive  index  for  a  frequency  change  of 
42%.  The  ease  and  low  cost  of  fabrication  of  the  metallic  PC 
vis-il-vis  a  dielectric  PC  and  metamaterials  make  them  ideal 
for  a  wide  range  of  applications.  Precise  control  over  the 
geometry,  the  choice  of  mode,  and  the  scalability  to  submi¬ 
crometer  dimensions  of  PCs  shows  promise  for  applications 
from  microwave  to  optical  frequencies. 
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lethal  challenge.  Although  the  toxicity  of 
MDMA  in  mice  is  less  consistent  than  in 
other  species7,  a  25%  lethal  dose  (LD25)  is 
estimated  as  about  50  mg  of  drug  per  kg  body 
weight8.  Consistent  with  this,  wild-type  ani¬ 
mals  given  increasing  doses  of  MDMA  under 
defined  conditions  succumbed  within  4  h  of 
drug  administration,  with  one  in  five  of  the 
animals  dying  at  a  dose  of  20  mg  kg-1  and 
three  in  ten  dying  at  40  mg  kg”  (Fig.  2c). 
However,  none  of  the  UCP-3_/_  animals 
died  after  administration  of  either  dose 
(n  =  5  and  n  =  1 0,  respectively) . 

Our  results  indicate  that  UCP-3  is  required 
for  the  rise  in  skeletal  and  core  temperature 
that  is  associated  with  MDMA  administra¬ 
tion.  Endogenously  produced  reactive  alde¬ 
hydes  can  also  stimulate  UCP-3  activity9,  and 
there  are  some  structural  similarities  between 
this  class  of  molecule  and  MDMA  and/or  its 
metabolites7.  Further  investigation  is  needed 
to  determine  whether  ecstasy  can  directly 
stimulate  uncoupling  activity. 

The  temperature  response  in  individuals 
who  take  excessive  amounts  of  MDMA 
varies,  and  this  variation  could  relate  to  the 
uncoupling  activity  of  their  skeletal  muscle. 
For  example,  it  is  known  that  thyroid-hor¬ 
mone  treatment  can  exacerbate  the  hyper¬ 
thermic  effects  of  amphetamines10  and  that 
UCP-3  is  upregulated  by  thyroid  hormone11. 
It  is  possible  that  UCP-2  and/or  UCP-3  could 
also  mediate  the  toxicity  of  agents  such  as 
ephedrine,  methamphetamine  and  cocaine, 
recreational  stimulants  that  can  produce  a 
similar  hyperthermic  response  and  which 
are  related  to  MDMA.  Our  demonstration 
that  UCP-3  is  a  molecular  mediator  of  the 
thermogenic  response  to  ecstasy  suggests 
that  agents  designed  to  target  uncoupling 
proteins  could  provide  the  basis  for  an 
important  new  therapeutic  direction. 
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Photonic  crystals 

Imaging  by  flat  lens 
using  negative  refraction 

The  positive  refractive  index  of  conven¬ 
tional  optical  lenses  means  that  they 
need  curved  surfaces  to  form  an  image, 
whereas  a  negative  index  of  refraction  allows 
a  flat  slab  of  a  material  to  behave  as  a  lens  and 
focus  electromagnetic  waves  to  produce  a  real 
image1.  Here  we  demonstrate  this  unique 
feature  of  imaging  by  a  flat  lens,  using  the 
phenomenon  of  negative  refraction  in  a  pho¬ 
tonic  crystalline  material.  The  key  advance 
that  enabled  us  to  make  this  observation  lies 
in  the  design  of  a  photonic  crystal2,3  with  suit¬ 
able  dispersion  characteristics  to  achieve  neg¬ 
ative  refraction  over  a  wide  range  of  angles. 

Although  negative  refraction  has  been 
demonstrated  at  microwave  frequencies  in  a 
quasi-homogeneous  metamaterial4,  imaging 
by  a  flat  lens  is  severely  constrained  because 
of  large  dissipation5,6  and  anisotropy  in  the 
metamaterial.  Plane-wave  negative  refraction 
at  specific  incident  angles  has  been  demon¬ 
strated  at  microwave  frequencies  by  using  a 
metallic  photonic  crystal  prism7  and  a  dielec¬ 
tric  photonic  crystal8.  However,  to  focus  a 
diverging  beam  from  a  point  source,  the 
material  must  exhibit  all-angle  negative 
refraction2  as  well  as  low  absorption. 

Figure  la  shows  the  image  of  a  microwave 
point  source  of  frequency  9.3  GHz  (wave¬ 
length,  3.22  cm)  placed  2.25  cm  from  a 
two-dimensional  flat  lens  made  of  a  photonic 
crystal  fabricated  from  an  array  of  cylindrical 
alumina  rods  (see  supplementary  informa¬ 
tion)  .  On  the  far  side,  a  high-quality  image  is 
seen  at  a  distance  of  2.75  cm.  Note  that  there  is 
an  image  of  similar  size  for  the  sub- wavelength 
source.  This  image  is  observed  only  in  a  narrow 
frequency  range,  between  9.0  and  9-4  GHz, 
with  the  best  focus  at  9.3  GHz.  Outside  this 
narrow  band,  at  all  other  frequencies  between 
2  and  1 2  GHz,  a  single  focus  point  is  not  seen. 

These  observations  can  be  understood  by 
examining  the  band  structure  of  the  photonic 
crystal  and  the  corresponding  equi-frequency 
surfaces,  from  which  an  effective  refractive 
index,  n^,  can  be  defined  and  calculated9.  In 
our  geometry,  the  central  axis  of  the  diverging 
beam  is  along  the  (1,0)  direction  of  the  square 
lattice  crystal.  Just  below  the  top  of  the  second 
band,  between  9.0  and  9.4  GHz,  the  equi-fre¬ 
quency  surfaces  move  inwards  with  increasing 
frequency,  consistent  with  negative  n,.n,  The 
condition  for  negative  refraction,  ri^fO)  <  0, 
holds  for  a  diverging  beam  of  sufficiently  large 
angle,  enabling  the  image  to  be  formed.  Figure  1 
therefore  also  demonstrates  wide-angle  nega¬ 
tive  refraction  by  a  photonic  crystal. 

The  value  of  nI:fris  necessarily  angle-depen- 
dent.  Inside  the  crystalline  lens,  the  electro¬ 
magnetic  field  is  highly  modulated  and 
simple  ray  diagrams  applicable  in  homo¬ 
geneous  media  cannot  be  used.  However, upon 


.  Figure  1 1maging  by  a  flat  lens.  a.  Microwave  efectric-field-intensity 
map  on  a  cross-sectional  view  of  the  two-dimensional  source-image 
(left-right)  system  for  the  flat  lens,  b,  Displacement  of  the  image  by 
4  cm  when  the  source  is  moved  up  by  the  same  distance.  Both 
panels  correspond  to  dimensions  of  37.5  x  30.0  cm.  The  intensity 
scale  on  the  source  side  varies  from  -  20  dB  (yellow)  to  -  48  dB 
(black),  and  on  the  image  side  from  -  30  dB  to  -  75  dB. 

emerging,  the  waves  create  a  focus  of  the  diver¬ 
gent  beam  emanating  from  the  point  object. 

Conventional  optical  systems  have  a  single 
optical  axis  and  limited  aperture,  and  cannot 
focus  light  onto  an  area  smaller  than  a  square 
wavelength1 .  In  contrast,  our  flat  lens  does  not 
have  a  unique  optical  axis  and  is  not  restricted 
by  aperture  size.  We  have  demonstrated  both 
of  these  features  by  moving  the  source  up  by 
4  cm:  the  image  moves  a  corresponding 
distance  in  the  same  direction  (Fig.  lb) . 

The  unique  properties  of  our  flat  lens 
provide  new  perspectives  on  imaging.  A 
particular  advantage  of  the  photonic  crys¬ 
talline  material  is  its  scalability  to  submicro¬ 
metre  dimensions  for  possible  applications 
from  microwave  to  optical  frequencies. 
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From  the  example  in  the  previous  section,  enhancing  efficiency 
for  convolution  in  the  presence  of  nonlinear  terminations  is  crucial, 
since  it  takes  up  most  of  the  run  time.  Although  we  only  confine 
ourselves  to  the  S-parameter  application  here,  other  parameters 
such  as  the  Z  or  Y  matrix,  preferred  for  RLC  circuits,  can  also  be 
used  easily.  Moreover,  this  algorithm  can  be  directly  incorporated 
into  multiconductor  transmission  lines  (MTL)  terminated  with 
nonlinear  loads  [1]  as  well. 
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ABSTRACT:  The  imaging  of  the  near  field  by  a  slab  of  negative  permit¬ 
tivity  medium  is  discussed.  The  performance  of  this  fiat  lens  is  limited  by 
the  working  wavenumber  k„  and  the  slab  thickness  d.  For  k,/l  <SC  0.937, 
the  resolution  of  the  flat  lens  can  be  smaller  than  the  wavelength  and  de¬ 
pends  logarithmically  on  /A,/l.  The  electric  field  is  finite  everywhere,  but 
will  diverge  algebraically  as  one  reduces  kod  -»  0.  ©  2003  Wiley  Periodi¬ 
cals,  Inc.  Microwave  Opt  Technol  Lett  39:  282-286,  2003;  Published  on¬ 
line  in  Wiley  InterScience  (www.interscience.wiley.com).  DOI  10.1002/ 
mop.  11 191 

Key  words:  near-field  imaging;  evanescent  waves;  negative  permittivity; 
refraction;  transmission 

1.  INTRODUCTION 

The  idea  of  negative  refraction  came  to  life  when  Veselago  [1] 
pointed  out  the  possible  existence  of  negative-index  media  (NIM) 
without  violating  Maxwell’s  equations.  Based  on  the  idea  of  neg¬ 
ative  refraction,  Pendry  [2]  proposed  that  a  flat  slab  of  such  a 
medium  can  focus  an  object  on  the  other  side  of  the  slab.  Not  only 
the  far  field  (which  constitutes  of  propagating  plane  waves),  but 
also  the  near  field  of  an  object  can  be  focused  through  such  a  slab. 
This  enables  one  to  circumvent  the  diffraction  limited  resolution 
and  allows  for  subwavelength  imaging.  The  proposal  of  this  “per¬ 
fect  lens”  led  to  growing  interest  and  also  initiated  some  contro¬ 
versy  [3-5].  It  is  now  becoming  clear  that  the  idea  of  a  perfect  lens 

PACS  numbers:  78.20.Ci,  41.20.Jb,  42.25.-p,  42.30.-d. 


Figure  1  Sketch  of  a  flat  lens  with  e  =  —  1,  p,  =  1,  and  slab  thickness 
d 


is  physically  sound  [6,  7]  and  also  experimentally  feasible  in 
systems  such  as  metamaterials  [8,  9]  and  photonic  crystals  [10]. 

Pendry  also  noticed  that  a  flat  slab  of  material  which  has 
negative  permittivity  (e  <  0)  can  act  as  a  lens  to  focus  the  near 
field.  The  near-field  optical  technique  has  been  pushed  to  have  a 
spatial  resolution  of  about  20  nm  [1 1]  at  optical  frequencies.  In  this 
paper,  we  discuss  the  near-field  imaging  of  this  lens  with  negative 
permittivity  e  <  0.  To  have  the  best  focus  of  the  near  field,  it  is 
desirable  for  the  slab  to  have  permittivity  e(to)  =  - 1  at  a  certain 
frequency.  As  required  by  causality,  no  material  can  have  static 
permittivity  e(ro  =  0)  =  —  1.  Only  at  a  certain  nonzero  frequency 
oj0  and  nonzero  wave  number  k0  =  m0/c,  s  =  —  1  can  this  be 
realized.  This  will  happen  when  one  crosses  a  resonant  frequency 
in  the  case  of  plasma  e[tu)  =  1  -  co2/(o2,  with  aip  the  plasma 
frequency.  This  flat  lens  with  e  =  —  1  will  amplify  the  P-polarized 
evanescent  waves,  thus  the  diffraction  limit  of  the  ordinary  lens 
which  only  focuses  propagating  waves,  will  be  broken.  In  order  to 
focus  the  S-polarized  near  field,  one  needs  material  with  negative 
permeability  p(a>)  =  —1.  As  pointed  by  Landau  [12],  at  very  high 
frequencies  such  as  at  optical  frequencies,  p  — »  1.  So,  a  fiat  lens  of 
any  medium  is  not  likely  to  focus  the  S-polarized  near  field  in 
optical  frequencies. 

The  presence  of  nonzero  k0  will  put  strong  restrictions  on  the 
performance  of  the  lens.  As  k0  — >  0,  the  lens  will  focus  more 
details  of  the  near  field.  In  this  paper,  we  discuss  in  depth  the 
limitation  of  the  superlens  due  to  finite  k0.  Shen  and  Platzman 
considered  the  imaging  of  a  point  dipole  [13].  Here,  we  give  a 
detailed  discussion  through  the  example  of  the  imaging  of  two- 
spike  source,  which  was  also  discussed  by  Pendry  [2],  The  diver¬ 
gence  of  the  field  in  the  limit  k0  — »  0  will  also  be  discussed. 

2.  SLAB  REFRACTION  OF  P-POLARIZED  PLANE  WAVES 

To  begin  with,  we  consider  a  plane  wave  with  P-polarization.  The 
incident  electric  field  is  given  by 

E0(kx)  =  e‘k'k+lk*(kxz  -  k,t).  (1) 

Let  the  thickness  of  the  slab  be  d.  In  our  coordinate  system,  the  z 
axis  is  along  the  direction  from  a  vacuum  to  the  slab  and  the  x  axis 
is  along  the  interface,  as  shown  in  Figure  1.  We  place  a  source  at 
z  =  0,  the  first  interface  of  the  slab  at  z  =  d/2,  and  the  second  at 
z  =  3d/2.  At  this  moment,  we  treat  the  slab  as  ordinary  dielec¬ 
trics.  The  field  configuration  of  the  slab  refraction  can  be  obtained 
using  standard  methods  [14].  The  reflected  and  transmitted  electric 
fields  are 
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E r(kx)  =  /f /.e'*'I+ d (~kxz  -  k7x), 

E.fc)  =  TPeik’x+>u,-d){kJi.  -  kzx), 

with  the  reflection  and  transmission  coefficients 

_  (1  -  T^r  (1  -  T2)7’ 
Rp  i _ i _ Tifi  » 


and 


(2) 

(3) 

(4) 


8=k  / k„ 

x  o 


T  =  (ekz  -  kn)/{ekz  +  kn),T=  exp(tM)-  (5) 

Here  kz  =  \/k^  —  kx  is  the  z  component  of  the  wave  vector  in 
vacuum  and  krz  =  V n^k^  —  kx  is  that  inside  the  slab.  The  electric 
field  inside  the  slab  is 


eik,i+lkpll2+lk,i7.-iin) 

■ 


i  -r 


(1  +  Ttleilni~d~li))z 


-(1  +T)(1 

-T 


•  (6) 


The  propagating  plane  waves  will  become  evanescent  when 
kx  >  k0,  so  that  k7  is  purely  imaginary,  that  is,  kz  =  i  Vkx  -  k^. 
These  near-held  waves  cannot  propagate  and  will  decay  exponen¬ 
tially  around  the  object  in  the  vacuum  along  the  z  axis.  For  those 
evanescent  waves,  the  magnetic  field  is  exponentially  small,  com¬ 
pared  with  the  electric  field.  Thus,  practically,  there  is  only  electric 
field,  as  in  the  electrostatic  case.  The  energy  density  of  the  incident 
evanescent  plane  wave  around  the  source  in  the  vacuum  is  u  = 
e2'*!74ir. 

The  above  expressions  of  the  slab  refraction  is  valid  for  arbi¬ 
trary  e  and  p..  For  NIM  with  e  =  - 1  and  p,  =  —  1,  since  krz  =  kz, 
the  denominator  of  T  in  Eq.  (5)  is  zero;  thus,  T  The 

evanescent  wave  will  be  amplified  since  the  transmission  coeffi¬ 
cient  TP  =  exp (Vk2  -  kid)  will  grow  exponentially  with  the 
slab  thickness  and  increasing  kx.  This  amplifying  factor  will  com¬ 
pensate  the  decay  in  the  vacuum  such  that  at  a  distance  d/2  away 
from  the  second  surface  of  the  slab,  Tpexp{ikzd )  =  1 .  Hence  all 
the  evanescent  waves  will  be  amplified  through  the  slab  and  a 
perfect  image  of  the  object  at  z  =  0  is  obtained  at  z  =  2d. 


3.  PROPERTIES  OF  TRANSMISSION  AND  REFLECTION 
COEFFICIENTS  FOR  EVANESCENT  WAVES 

Here  we  consider  a  slab  with  e  =  —  1  and  ju,  =  1 .  A  slab  of  metal 
at  certain  frequencies  will  have  this  property.  Although  this  lens 
will  not  perfectly  focus  all  the  evanescent  waves,  it  can  recover 
part  of  the  evanescent  waves  under  certain  conditions.  The  slab 
will  still  be  able  to  achieve  subwavelength  imaging.  Since  in  this 
case  the  index  of  refraction  is  purely  imaginary  with  n  =  i,  the  z 
component  of  the  wave  vector  for  the  evanescent  wave  is  krz  = 
i  \Zk-i  +  kf,.  To  proceed,  we  define  S  =  kjk0,  which  is  greater 
than  unity  for  evanescent  waves.  From  the  expression  of  F  in  Eq. 
(5),  one  obtains  T  =  — (82  +  VS4  —  1)  <  —  1.  Unlike  in  the  case 
for  NIM  with  e  =  —  1  and  ju,  =  —  1,  here  T  is  finite. 

We  first  look  at  the  property  of  Tp.  From  the  expression  in  Eq. 
(3)  of  Et(kx),  one  can  see  that  if  Tpexp(ikzd)  =  I,  then  E t(kx) 
with  z  >  2d  will  be  exactly  the  same  as  E 0(kx)  with  z  a  0.  Here 
the  factor  exp (ikzd)  is  the  damping  factor,  due  to  the  propagation 
of  the  evanescent  waves  in  the  vacuum.  Since  r  <  0,  TP  will 
diverge  at  TT  =  —  1.  This  condition  can  be  re-expressed  in  terms 
of  8  in  the  following  equation: 


Figure  2  Plot  of  Tpexp(ikzd)  for  (a)  k0d  =  1.0,  (b)  k0d  =  0.1 


M  =  [ln(S2  +  VS4  -  1 )]/  VS2  +  1.  (7) 

The  maximum  of  the  function  on  the  right  side  of  above  equation 
is  about  0.937  261  at  S  “  2.356  574.  So  if  k0d  >  0.937,  - 1  < 
T  exp (ikrzd)  <  0,  there  is  no  singularity  in  Tp  with  Tpexp(ikzd) 
being  negative,  as  can  be  seen  in  Figure  2(a)  with  k0d  =  1 .0.  Even 
though  the  maximum  of  \Tpexp(ikzd)\  can  be  much  larger  than  1, 
\T pcxp(ik7d)\  is  not  very  flat.  Evanescent  waves  with  different  kx 
will  have  different  amplifying  factors.  For  a  source  with  many 
component  of  kx,  the  transmitted  evanescent  waves  will  be  se¬ 
verely  distorted.  So  there  is  no  subwavelength  imaging  for  k0d  > 
0.937.  On  the  other  hand,  if  k0d  <  0.937,  Eq.  (7)  will  give  two 
solutions,  8,  and  S2,  for  8  with  8[  <  82.  If  k0d  <§:  0.937,  one  has 
8X  — >•  1  and  S2  —  —2  In (k0d)/k0d.  If  kx  >  k™*  =  S2k0  -  -2 
In (k0d)/d,  except  in  the  vicinity  of  k‘"ax.  Tpexp(ikzd)  decays 
exponentially  to  zero  with  increasing  kx.  More  specifically,  for  kx 
»  k'"ax.  one  has  FT  — »  0,  so  that  Tp  ~  4 S4e~txd.  Details  of  the 
object  with  high  Fourier  components  k  >  k rx  will  be  lost  in  the 
imaging  process.  What  is  amplified  through  the  slab  is  a  finite 
range  of  Fourier  components,  kx  £  (&™in,  kx'ax)  with  k'x>n  =  S^q. 
Within  this  range,  except  in  the  vicinity  of  k'x'"  and  &"’0X, 
Tpexp(ikzd)  is  very  flat  and  Tpexp(ikzd)  —  1 .  This  can  be  seen  in 
Figure  2(b)  with  k0d  =  0.1.  In  the  neighborhood  of  k0d  ~  0.937, 
the  lens  cannot  focus,  even  though  some  components  of  the  eva¬ 
nescent  waves  do  become  amplified  through  the  slab. 

The  reflection  coefficient  RP  has  similar  behavior  as  TP.  From 
the  expression  of  TP  and  RP,  the  difference  between  them  is  T P  — 
RP  =  (T  -  r)/(l  -  r7j.  Since  TT  <  0,  there  is  no  singularity 
in  TP  -  RP.  The  difference  will  diverge  exponentially  for  eva¬ 
nescent  waves  with  kx  —*  <».  The  behavior  of  RPexp(ikzd)  for 
different  values  of  k0d  is  shown  in  Figure  3.  For  finite  k0d  <SC 
0.937,  RPcxp(ikzd)  is  flat  with  nonzero  values  for  £"’in  <5C  kx  <3C 
&™ax.  This  is  due  to  the  mismatch  of  the  surface  impedance 
between  the  vacuum  (Ve/g  =  1)  and  the  slab  (Vs/g  =  i). 
Unlike  the  behavior  of  TP ,  which  decays  exponentially  for  large 
kx,  the  amplitude  of  RP  increases  as  I/?  J  ~  282  for  kx  »  k""'x 
with  kgd  «  0.937. 
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4.  NEAR  FIELD  IMAGING  OF  A  TWO-SPIKE  OBJECT 

To  illustrate  the  performance  of  the  above  lens,  we  take  as  an 
example  the  two-spike  object  [2],  The  electric  field  is 


E, 


=  |  dkxvk,e'kzZ(& 
Jo 


cos  k^x  —  ix  —  sin  k,x 

kx 


with  the  coefficient  vk  given  by 

vlx  =  (4/irkj)sin  kxa  cos  kxb. 


(8) 


(9) 


Figure  4  The  near  field  of  the  two-spike  object,  E0z  (top)  and  E0x 
(bottom)  with  a  =  0.2,  b  =  1,  d  =  1,  and  k0d  =  0.1 


En  =  —  S)l  V(z-  ('x)arctan - r, — —rr  , 

7r  z  -  i(x  ±  b) 


E,  =  2  (2  -  «)arctan  jZ^JFTb)  ’ 


(10) 


Here,  2b  is  the  distance  between  the  two  spikes  and  2 a  is  the 
individual  spike  width.  The  waves  will  propagate  freely  in  the 
vacuum  for  kx  <  kQ.  Since  kz  is  purely  imaginary  for  evanescent 
waves  with  kx  >  k0,  the  near  field  part  of  E0  of  the  two-spike 
object  is  real  as  expected,  which  is  illustrated  in  Figure  4. 

The  transmitted  near  field  can  be  obtained  by  the  integration  of 
Eq.  (3)  with  vk  /kx  over  kx.  For  the  source  shown  in  Figure  4,  the 
transmitted  near  field  is  shown  in  Figure  5  with  k0d  =  0.1.  The 
reflected  field  and  the  field  inside  the  slab  can  also  be  obtained 
similarly.  The  near  field  distribution  in  the  vicinity  of  the  slab  with 
k0d  =  0. 1  is  shown  in  Figure  6.  One  notices  that  there  is  quite  a 
lot  of  field  accumulated  around  the  two  surfaces.  The  amplitude  of 
the  field  around  the  two  surfaces  reaches  about  40,  as  compared 
with  that  of  the  source,  which  reaches  1  (this  will  be  explained 
subsequently). 

The  imaging  performance  for  slabs  with  different  values  of  k0d 
is  shown  in  Figure  7.  One  can  see  that,  as  the  value  of  k0d  is 
increased,  the  two  spikes  at  the  imaging  plane  z  =  2d  will  be  less 
distinguishable.  When  the  value  of  k0d  is  close  to  and  beyond 
0.937,  the  slab  lens  loses  the  power  to  amplify  evanescent  waves 
and  the  two  spikes  disappear. 

In  order  to  address  the  issue  of  singularity  of  the  near-field 
imaging,  we  consider  the  electrostatic  limit  ka  — »  0.  In  this  limit, 
we  assume  that  e  =  -1  still  holds,  even  though  the  law  of 
causality  will  be  violated.  Subsequently,  one  has  k7  =  krz  =  ikx. 
As  a  result,  r  — >  °°,  and  we  obtain 


through  a  slab  of  e=-1 


X 


through  a  slab  of  e=-1 


x 


Figure  5  Transmitted  near  field  of  the  two-spike  source  (Fig.  4),  En 
(top)  and  Elx  (bottom),  for  z  &  2d 
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eE. 


Figure  6  Near-field  e£z  (top)  and  Ex  (bottom)  in  the  vicinity  of  the  slab 
for  the  two-spike  object  with  a  =  0.2,  b  =  1,  d  =  1,  and  k0d  =  0.1 


Here,  we  use  the  following  integral: 

f"  sin  ka  ,,  a 

I  dk — — e  *£  =  arctan^,  (11) 

*  o 

for  any  complex  f  with  >  0. 

The  above  expression  for  E,  is  only  valid  for  z  &  2d.  The 
transmitted  electric  field  E,  diverges  in  the  region  with  z  E  (3d/ 2, 
2d).  The  reflected  field  Er  diverges  for  z  E  (0,  d/2).  The  field  Es 
inside  the  slab  diverges  everywhere  except  on  the  line  z  =  d.  The 
field  Es  on  this  line  is 

Es  =  ^2>- <x)arctan  W_-(x±b)  • 

This  means  that  one  has  a  perfect  image  of  the  source  in  the  middle 
of  the  slab  besides  another  one  outside  the  slab.  This  scenario  is 
similar  to  the  case  of  far-field  imaging,  except  that  here  one  has 
subwavelength  imaging. 

For  any  physical  material,  the  dielectric  constant  should  always 
be  positive  as  one  approaches  the  electrostatic  limit.  The  presence 
of  nonzero  k0  for  e  =  —  1  will  remove  the  divergence  which  is  due 
to  an  accumulation  of  field  along  the  two  surfaces  of  the  slab.  For 
small  but  finite  k0d,  the  field  is  bounded  and  finite  everywhere.  As 
one  decreases  k0d  — >  0,  the  electric  field  will  diverge.  To  check 


the  divergence,  we  consider  for  example  the  transmitted  field  Er  at 
z  =  d/2,  where  the  field  is  more  likely  to  diverge.  Since 

f-  2lnM  g3f/2  A o 

|Er|~  ~  (MT3, 

I  •/  -2  In  tod 

the  electric  field  around  the  first  surface  of  the  slab  will  diverge  as 
E  ~  (kod)~3.  For  the  field  around  the  second  surface,  we  consider 
Et  at  z  =  3d/2.  One  has 


The  electric  field  around  the  second  surface  of  the  slab  will  diverge 
as  E  ~  ( k0d )_1. 

5.  CONCLUSION  AND  DISCUSSION 

In  summary,  we  discussed  the  limitation  of  the  lens  made  of  a  slab 
with  negative  permittivity.  For  small  k^d,  the  lens  can  focus 
partially  subwavelength  details  of  the  P-polarized  near  field  of  an 
object.  The  resolution  of  the  slab  lens  depends  logarithmically  on 
(k0d)~  \  since  k'”“x  —  -~2  In (k0d)ld.  For  fixed  k0,  in  order  for  the 
lens  to  pick  up  the  high  Fourier  component  of  the  evanescent 
waves,  one  may  reduce  the  slab  thickness  d  such  that  k0d  <K 
0.937.  Given  the  fact  that  the  minimum  thickness  of  a  material 
under  normal  condition  is  about  d  ~  10“ 8  cm,  the  working 
frequency  of  the  slab  lens  will  be  ai/2ir  <  5  X  1017  Hz,  which 
includes  the  optical  frequency  range.  At  w/2ir  =  3  X  1014  Hz  with 
slab  thickness  d  =  100  nm,  one  has  k0d  =  0.1  and  —  10s 
m-1,  and  a  spatial  resolution  of  about  60  nm  can  be  obtained.  In 
the  limit  k0d  — >  0,  all  the  evanescent  waves  will  be  amplified; 
thus,  the  slab  lens  will  become  perfect  for  near-field  imaging.  The 
electromagnetic  field  is  finite  everywhere  but  will  diverge  alge¬ 
braically  in  the  limit  k0d  — >  0  on  the  surfaces  of  the  slab  lens. 
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Figure  7  Images  of  evanescent  waves  in  the  plane  at  z  =  2d  with  a  = 
0.2,  b  =  I ,  and  d  =  1  for  different  values  of  k0d  =  10-4,  10-3,  10-2, 
and  10_1 
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ABSTRACT:  A  high-power  dual-wavelength  laser  is  constructed  in  two 
fiber-loop  mirror  cavities  by  utilizing  a  common-gain  medium.  The  output 
power  and  side-mode  suppression  ratio  (SMSR)  of  the  lasers  were  investi¬ 
gated  for  a  bidirectionally  pumped  system.  The  outcome  of  the  experiment 
shows  that  two  lasing  signals  can  be  obtained  at  the  wavelengths  of 
1546.04  and  1550.64  nm  at  a  less  than  1-dB  peak  power  difference,  at  peak 
powers  of  7.30  and  6.57  dBm,  respectively.  This  result  was  obtained  with¬ 
out  any  intracavity  attenuators.  The  side-mode  suppression  ratio  of  both 
signals  is  around  46  dB.  ©  2003  Wiley  Periodicals,  Inc.  Microwave  Opt 
Technol  Lett  39:  286-287,  2003;  Published  online  in  Wiley  InterScience 
(www.interscience.wiley.com).  DOI  10. 1002/mop.  11 192 
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Figure  1  Configuration  of  the  dual  wavelength  bidirectional  pumped 
dumbbell  EDFL 

INTRODUCTION 

Multiple-wavelength  optical  sources  have  been  investigated  as  one  of 
the  promising  technologies  for  applications  in  wavelength  division 
multiplexed  (WDM)  systems.  Many  techniques  have  been  employed 
in  designing  a  laser  system  emitting  multiple  wavelengths  using 
erbium-doped  fiber  laser  (EDFL).  Park  et  al.  [1]  used  eight  fiber-loop 
mirrors  (FLMs)  to  produce  eight  wavelengths.  However,  the  output 
power  of  all  the  channels  is  not  equal  and  is  relatively  low.  Another 
technique  to  produce  multiple  wavelengths  involves  the  use  of  a  twin 
core  erbium-doped  fiber  (EDF)  [2].  The  output  wavelengths  are 
selected  by  using  fiber  Bragg  gratings  (FBGs).  However,  the  EDF 
used  is  a  highly  specialized  fiber  that  is  costly,  and  the  output  power 
reported  with  this  technique  is  relatively  low.  Another  approach  is  to 
use  a  periodic  filter  to  produce  multiple  wavelengths  [3],  To  prevent 
single-frequency  lasing,  a  frequency  shifter  is  introduced  in  the  ring 
cavity.  It  requires  two  gain  media  to  compensate  the  loss  introduced 
by  the  frequency  shifter.  This  approach  is  rather  complicated.  A 
high-power  output  dual- wavelength  laser  is  also  reported  in  [4], 
However,  the  output  is  taken  from  two  fibers,  instead  of  one.  Re¬ 
cently,  a  linear  cavity  configuration  using  cascaded  FBG  and  EDF  has 
been  reported  [5],  The  output  power  of  the  multiple  lasers  is  consid¬ 
erably  high;  however,  this  technique  does  not  deploy  a  single  EDF 
coil.  Furthermore,  the  reflectivity  of  the  FBGs  are  used  to  create  an 
equal  net  gain  of  the  laser  system  at  multiple  wavelengths. 

In  this  paper,  a  new  and  simple  configuration  utilizing  FLMs  with 
an  intraloop  tunable  bandpass  filter  (TBF)  is  proposed.  The  signals 
share  most  of  the  cavities,  thus  reducing  the  cavity  discrepancy.  This 
includes  a  common  section  of  EDF  to  provide  amplification  for  both 
lasing  wavelengths.  In  this  configuration,  the  two  signals  have  almost 
equivalent  characteristics  in  terms  of  power  and  side  mode  suppres¬ 
sion  ratio  (SMSR)  at  certain  wavelengths.  The  1-dB  difference  peak 
power  produced  by  this  system  is  higher  than  previously  reported 
results  for  a  system  employing  relatively  low  pump  power,  and  short 
EDF  with  low-ion  concentration. 

EXPERIMENTAL  CONFIGURATION 

The  setup  for  this  experiment  follows  a  dumbbell-shaped  config¬ 
uration,  as  shown  in  Figure  1.  The  left  side  of  the  dumbbell  has 
been  modified,  such  that  it  consists  of  two  FLMs.  The  amplified 
spontaneous  emission  (ASE)  signal  will  only  pass  through  it  in  one 
direction  via  a  three-port  circulator.  The  directivity  loss  for  port 
3-1  is  around  65  dB,  which  is  enough  to  prevent  the  signal  from 
taking  the  shortcut  through  it.  Two  other  ports  have  1-dB  insertion 
loss  each.  Therefore,  the  ASE  signal  enters  the  loop  from  port  2-3 
and  propagates  through  the  TBF  and  out  of  the  loop  via  port  1-2. 
A  fused  fiber  coupler  is  used  to  obtain  the  laser  output  at  the  20% 
port.  The  EDF  used  is  10-m  long  with  a  core  diameter  of  2  ptm  and 
erbium  concentration  of  440  ppm.  Both  of  the  TBFs  have  a  3-dB 
passband  of  1  nm  and  tunable  from  1525  nm  to  1565  nm. 

On  first  pass,  the  ASE  is  equally  split  into  the  two  loops.  In 
each  loop,  most  of  the  ASE  signal  will  be  filtered  out,  and  only  the 
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I.  INTRODUCTION 

The  existence  of  a  medium  with  a  negative  (n<0)  index 
of  refraction,  raised  several  years  ago  [1],  has  recently  been 
demonstrated  experimentally  [2].  One  of  the  most  striking 
properties  of  negative  index  materials  (NIMs)  is  the  negative 
refraction  for  plane  waves  across  the  interface  between 
positive-index  material  (PIM)  and  a  NIM.  Negative  refrac¬ 
tion  means  that  when  radiation  passes  through  an  interface 
between  a  PIM  and  a  NIM,  the  refracted  beam  is  on  the  same 
side  of  the  normal  as  the  incident  beam  (see  Fig.  1),  in  con¬ 
trast  to  the  usual  positive  refraction  in  which  they  are  on 
opposite  sides  of  the  normal. 

In  studies  of  negative  refraction,  it  is  essential  to  represent 
incident  waves  as  localized  wave  packets,  rather  than  plane 
waves,  since  all  physical  sources  of  electromagnetic  waves 
produce  radiation  fields  of  finite  spatial  and  temporal  extent 
because  the  sources  are  always  of  finite  spatial  extent  and 
because  they  only  radiate  for  a  finite  time.  Hence  treatments 
of  this  problem  which  study  waves  that  extend  over  infinite 
distance  in  all  or  some  directions  cannot  be  trusted  to  reli¬ 
ably  predict  the  direction  in  which  a  wave  will  be  refracted, 
and  in  fact  treatments  based  on  such  extended  waves  [3] 
have  led  to  a  direction  of  refraction  opposite  to  that  which 
one  finds  for  spatially  localized  wave  packets,  resulting  in  a 
great  deal  of  controversy  and  confusion.  Although  several 
treatments  using  waves  of  infinite  extent  in  some  direction 
(e.g.,  a  plane- wave  front  [4])  have  obtained  negative  refrac¬ 
tion,  since  such  a  model  is  unphysical,  for  the  reasons  given 
above,  we  cannot  have  confidence  in  conclusions  obtained 
from  it. 

In  this  paper,  we  treat  refraction  of  a  localized  wave 
packet  at  a  PIM-NIM  interface  both  analytically  and  by 
simulations,  demonstrating  that  it  refracts  negatively.  We 
also  present  both  analytic  and  numerical  studies  of  wave 
packets  constructed  from  a  small  number  of  plane  waves. 
Our  purpose  in  doing  this  is  to  give  a  plausible  explanation 
for  why  the  two-plane-wave  model  studied  by  Valanju  et  al. 
[3]  gives  a  misleading  answer.  We  find  that  in  all  cases, 
including  the  model  of  Valanju  etal.,  the  energy  and  mo¬ 
mentum  of  the  wave  refract  negatively.  Since  electromag¬ 
netic  waves  are  detected  only  when  they  either  give  up  en¬ 
ergy  to  or  exert  a  force  on  a  detector,  the  relevant  direction  of 
propagation  to  consider  is  that  of  the  region  of  space  in 
which  the  energy  and  momentum  of  the  wave  are  nonzero. 


Without  sources,  Maxwell’s  equations  are  VD  0,  V 
XH  <9,D,  VXE+<9,B  0,  V  B  0.  For  plane  waves  of 
wave  vector  k  and  frequency  «,  only  three  equations  are 
independent.  Using  the  usual  relationships  between  D(r)  and 
E(r)  and  between  B(/)  and  H(f)  [5]  one  obtains  for  such 
plane  waves  kXH  we(w)E,  kXE  m/z(o>)H.  Combin¬ 
ing  these  equations  gives  us  a  functional  relationship  be¬ 
tween  to  and  k.  Wave  propagation  is  only  permitted  for 
( e,fi,n>0 )  or  (e,jz,n< 0)  [1],  In  the  latter  case,  (E,H,k) 
will  form  a  left-handed  triplet  while  in  the  former  case,  for 
an  ordinary  material,  (E,H,k)  will  form  a  right-handed 
triplet. 

In  Sec.  II,  we  treat  the  negative  refraction  of  wave  packets 
and  beams.  The  analysis  of  refraction  of  finite  number  of 
plane  waves  will  be  done  in  Sec.  in. 

II.  NEGATIVE  REFRACTION  OF  WAVE  PACKETS 

A  wave  packet  localized  in  a  compact  region  of  space,  as 
occurs  in  all  experimental  situations,  can  be  constructed  from 
a  continuous  distribution  of  wave  vectors.  Consider  such  a 
wave  packet  incident  from  outside  the  NIM,  E 
yE0jd2kf(k  k0)eWkr  "(k),)1  with  w(k)  ck.  Here  we 
only  consider  S-polarized  waves.  The  P-polarized  waves  can 
be  treated  similarly,  however.  Throughout  the  paper,  we 
choose  the  z  axis  from  PIM  to  NIM  normal  to  the  interface 
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FIG.  1.  Time-lapse  snapshots  of  the  electric-field  intensity  of  a 
propagating  Gaussian  wave  packet  refracting  negatively  at  a  PIM- 
NIM  interface.  Arrows  indicate  the  directions  of  motion.  The  center 
wave  number  is  k0  1/5  with  incident  angle  ir/6.  The  spatial  extent 
of  the  incident  wave  packet  is  Ax  A z  10.  The  time  step  is  50 
with  speed  of  light  c  1 .  The  dispersion,  Eq.  (8),  was  used  for  NIM 
and  n  p.  1  for  PIM. 
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and  the  x  axis  along  the  interface.  If  /( k  ko)  drops  off 
rapidly  as  k  moves  away  from  ko,  w(k)  can  be  expanded  in 
a  Taylor  series  to  first  order  in  k  k0  to  a  good  approxima¬ 
tion.  This  gives 

E  y£,0et'(k°'r  w(ko>'>]g(r  ygt),  (I) 

with  g(R)  Jd2kf(k  ko)e,(k  k°)  R  and  \g 
Vk&>(k)|k  kfl. 

Inside  the  NIM,  k  and  ko  in  the  argument  of  the  exponent 
get  replaced  by  kr  and  kr0  which  are  related  to  k  and  k0  by 
Snell’s  law, 

krx  kx,  krz  yj{nr<i>lc)2  k2x.  (2) 

Here  nr  is  the  refractive  index  for  the  NIM  and  is  a  function 
of  cd.  Then  the  wave  packet  once  it  enters  the  NIM  is  given 
by 

Er  y£0e[,(k'"'r  “<k'0«]gr(r  vr),  (3) 

where  gr(R)  fd2kf( k  k0)fke,R(kr  kro)  and  fk  is  the 
transmission  amplitude  for  an  incident  plane  wave  of  wave 
vector  k.  It  is  the  standard  expression  for  this  quantity  for  the 
two  polarizations  of  the  incident  plane  wave  [6].  Here  kr0 
denotes  kr  evaluated  at  k  ko  and  \gr  Vka>(kr)  evaluated 
at  kr  kr0.  Let  us  expand  kr  kr0  in  the  exponential  function 
in  the  expression  for  gr(R)  in  a  Taylor  series  in  k  k0  to  first 
order,  kr  kr0~(k  ko)-Vk(kr  kr0)|k  k().  Substituting 
this  in  the  expression  for  gr( r),  we  obtain  gr(R) 
f  d2kf(k  k0)rke'R  [(k  ko>-vA  k'»> I.  if  the  width  of  the 
distribution  of  wave  vectors  /(k  ko)  is  small  compared  to 
the  range  of  k  over  which  tk  varies  significantly,  we  can  to  a 
good  approximation  simply  evaluate  this  quantity  at  k  ko 
and  put  it  outside  the  integral  over  k.  Then,  the  transmission 
coefficient  of  the  wave  packet  is  simply  given  by  |rkJ2. 

If  we  carry  out  the  expansion  of  w(kr)  to  second  order  in 
k  ko,  we  are  able  to  show  that  the  wave  packet  spreads  out, 
but  if  the  length  and  width  of  the  packet  are  much  larger  than 
the  wavelength  corresponding  to  the  wave  vector  k0  at  the 
peak  in  /( k  k0),  we  find  that  the  amount  that  the  packet 
spreads  out  in  a  given  time  interval  is  much  smaller  than  the 
distance  traveled  by  the  packet  in  that  time.  Then  clearly 
under  such  reasonable  conditions,  the  wave  packet  will  re¬ 
main  sufficiently  well  defined  to  be  able  to  observe  the  re¬ 
fraction  of  the  packet.  The  expansion  of  the  frequency  in  a 
Taylor  series  is  valid  for  a  sufficiently  narrow  distribution, 
/(k  ko)- 

In  order  to  get  an  explicit  expression  for  g(R),  let  the 
wave  packet  have  a  Gaussian  form  /( k)  (AxAz/ v)e.xp 
[  a£(A xf  k2(Az)2].  Expanding  k,  in  a  Taylor  series  around 
kr0,  we  get 

*r(R)  exp[  C2/4(A*)2  C2/4(Az)2]  (4) 

with  Cx  Rx+(cnrlvr  l)(kx0/krz0)Rz,  Cz  ( cnr / 
vr)(kzo/krZQ)Rz,  and  vr  c(dnrcD/dcD)  '.  From  the  above 


expressions,  one  can  see  that  the  Gaussian  wave  packet 
moves  with  \gr .  Due  to  the  dispersion,  the  wave  packet  is 
deformed  in  the  NIM. 

A  NIM  is  dispersive  and  causality  demands  that 
d(B(o)/dca>  1  and  d(/icD)/da>>  1  for  nearly  transparent  me¬ 
dia  [5,7].  For  an  isotropic  NIM,  since  nr  is  a  function  of  cd 
only,  \gr  c(dnrcD/do)  ’(ckr/nrw)  vrkr  with  kr  the 
unit  vector  in  the  direction  of  kr .  Since  vr  is  always  positive 
for  transparent  media  as  required  by  causality,  the  group  ve¬ 
locity  will  be  refracted  opposite  to  the  direction  of  wave 
vector  kr . 

The  magnetic  field  obtained  from  the  electrical  field 
through  H  (l/w/u,)kXE  is 

Hr  -f  d2kf( k  k0)— -T7T(krXy)e'k'''r  -<*>'  (5) 

C  fXr\Kr) 

with  Eq  tk(E0,  from  which  we  find  the  Poynting  vector  to 
be 


Sr  Re  ErXRe  Hr 


d2k 


d2k'f( k  ko)/(k'  k0) 


nr(k) 

fir(k) 


Xcos[kr-r— w(^)r]cos[k'-r  w(ife')f]kr,  (6) 


where  we  have  used  the  fact  that  kr-  y  0.  While  there  is  no 
question  that  the  Poynting  vector  at  a  point  in  a  medium 
gives  the  local  direction  of  energy  flow,  it  does  not  give  us 
the  direction  of  energy  flow  by  a  wave  packet  or  a  group  of 
plane  waves  as  a  whole  since  the  direction  of  the  Poynting 
vector  varies  with  space.  The  integral  of  the  Poynting  vector 
over  all  space,  Pr  / (Sr)dr,  however,  gives  the  total  mo¬ 
mentum  carried  by  a  wave  packet.  This  quantity  divided  by 
the  volume  over  which  the  wave  packet  is  nonzero  is  the 
average  of  the  Poynting  vector  over  the  whole  wave  packet. 
Either  way,  this  integral  clearly  represents  the  direction  of 
motion  of  the  wave  packet  in  the  medium.  From  the  above 
expression  for  Sr ,  one  has 


(|£'|2/2c) 


d2kf{  k  k0)2 


nr(k) 
Hr(k)  r 


(7) 


Let  us  consider  a  coordinate  system  whose  z  axis  is  along  k0. 
The  function  /( k  k0)2  will  then  be  a  function  of  kx  and  kz 
symmetrically  peaked  around  kx  0  and  kz  k0.  Then  writ¬ 
ing  Eq.  (7)  as 


P  r 


(K\2/2 c) 


d2kf( k  k0): 


nr(k)  kxx+knz 
fir{k)  |fcr| 


we  can  see  that  since  the  wave  number  k  is  an  even  function 
of  kx ,  the  integrand  is  an  odd  function  of  kx  and  hence  the  x 
component  vanishes.  Therefore,  Pr ,  which  as  argued  above 
represents  the  propagation  direction  of  the  wave  packet,  is 
opposite  in  direction  to  kr0  krz 0z,  i.e.,  in  the  direction  of 
the  group  velocity.  Hence,  the  energy  refracts  negatively. 
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FIG.  2.  Intensity  of  electric  field  Re  E  of  a  beam  with  kn  y[5 
and  the  Gaussian  weight  f(kx)  e  (lw±)  .  The  incident  angle  of 
the  beam  is  6  tt/12. 


The  negative  refraction  of  the  wave  packet  is  illustrated 
by  numerical  simulation  in  Fig.  1.  We  use  the  following  dis¬ 
persion  relation 

nr(co)  (l/w)V(w2  (»l)((o2  «2)/(w2  Wq)  (8) 

for  the  NIM  with  w0<  w<  cob .  The  permeability  is  fir 
(w2  (o2b)/(o)2  wq).  The  numbers  we  used  in  the  calcu¬ 
lation  are  <w0  1 ,  a>b  3,  o)p  \j\0,  and  c  1.  Figure  1 
shows  stroboscopic  snapshots  of  the  electric-field  intensity  of 
a  propagating  wave  packet  incident  on  a  PIM-NIM  interface 
[8].  The  negative  refraction  of  the  wave  packet  is  clearly 
evident. 

For  completeness,  let  us  consider  a  beam  given  by 

E  E0  zffc1ei(ko+k^)r/(fc±).  (9) 

Here  k±  is  perpendicular  to  k0  and  f(kL)  assumes  a  Gauss¬ 
ian  form.  Note  that  this  construction  is  different  from  that  of 
Kong  et  al.  [9]  and  Smith  et  al.  [10]  in  that  the  width  of  the 
incident  packet  is  made  finite  in  directions  perpendicular  to 
the  direction  of  propagation.  The  electric  field  E  of  the  beam 
is  shown  in  Fig.  2  [8].  Because  the  NIM  is  highly  dispersive, 
the  incident  beam  once  it  enters  the  NIM  will  no  longer  be  a 
beam.  It  will  be  a  localized  wave  packet  instead,  although  it 
is  difficult  to  see  this  in  the  figure.  Just  as  for  the  wave 
packet,  the  beam  intensity  also  refracts  negatively. 

III.  NEGATIVE  REFRACTION  OF  PACKETS 
CONSTRUCTED  FROM  A  FINITE  NUMBER 
OF  PLANE  WAVES 

Although  we  have  already  considered  the  refraction  of  a 
wave  packet  when  it  enters  a  NIM  from  a  PIM,  we  next 
consider  the  refraction  of  wave  packets  made  up  of  a  finite 
number  of  plane  waves.  Our  reason  for  doing  this  is  to  pro¬ 
vide  a  plausible  explanation  for  why  two  plane  waves  ex¬ 
ample  of  Valanju  et  al.  appears  to  give  positive  refraction. 


For  the  cases  of  two,  and  three  plane  waves  analytical  ex¬ 
pressions  are  obtained  for  the  Poynting  vector,  momentum, 
and  velocity  of  interference  pattern.  First  consider  the  case  of 
two  plane  waves  in  the  xz  plane  incident  from  PIM  to  NIM 
where  the  interface  is  at  z  0.  Let  wave  vectors  and  frequen¬ 
cies  be  (ki  ,o)[)  and  (k2 , oj2) -  We  set  the  polarization  in  the 
y  direction  as  before.  Suppose  Aw  w2  w, >0.  The  inci¬ 
dent  wave  in  PIM  is 


E  2 E0eiK  r  iS1'  cos(Ak-r/2  Awf/2)y  (10) 

with  K  (k)  +  k2)/2,  H  (w!  +  w2)/2,  Ak  k2  k^andEo 
the  wave  amplitude  of  each  plane  wave.  The  electric  field  of 
the  refracted  waves  is 

Er  2E'neiK'  r  iSlt  cos(Akr-r/2  Awf/2)y  (11) 


with  Kr  (krl  +  kr2)/2  and  Akr  k,2  krl ,  where  E'0 
tKE0 ,  krl  and  kr2  are  related  to  k,  and  k2,  respectively,  by 
Eq.  (2). 

The  relatively  long-wavelength  cosine  function  in  Eq. 
(11)  moves  in  the  NIM  with  a  velocity  [11] 

vr  (Aw/|Akr|2)(AkJcx+  Akrzz),  (12) 

assuming  that  |Ak|^|K|.  From  the  above  expression,  it  is 
evident  that  vrx>0  if  Akx>0.  Since  wi<w2,  we  have  0 
<«(w!)^n(w2)  and  nr(w[)<nr(w2)<0  by  the  require¬ 
ment  of  causality  which  requires  d(nr<o)/d(o>0.  One  has 

knz  k2riz  kn  k2r2  +  k2x  *i*>0-  Since  krz  |kj,  vrz 
>0,  the  group  refraction  appears  to  be  positive  [3].  This  is 
due  to  the  simple  fact  that  u„>0  if  vx>0.  Proper  dispersion 
must  give  vrz>0  since  the  energy  should  propagate  away 
from  the  interface.  But  we  shall  see  that  the  above  picture  is 
not  true  for  the  energy  flow. 

Let  us  determine  the  time-average  Poynting  vector  (S,.). 
Using  the  magnetic  field  corresponding  to  Er  of  Eq.  (11), 
Hr  EgS2  i(krjzx  kjXz)e‘kri'r  '“Wwj ,  (Sr)  is  defined  to 
be  the  time  average  of  Sr  over  the  period  corresponding  to 
the  average  frequency  of  the  two  plane  waves.  It  is  given  by 


1 

~ ( 1 +  cos 


2 

A0r)|Ei|2S 

j  1 


&  +  £&  ,  (13) 

W:  W: 


where  A</>r  Akr-r  A. cot.  Since  krz<0,  one  has  Sx<0  and 
Sz> 0.  Thus,  contrary  to  the  refraction  of  the  cosine  function 
in  Eq.  (11),  the  Poynting  vector  is  directed  in  the  negative 
refraction  direction,  i.e.,  refracts  negatively. 

We  shall  now  demonstrate  that  by  including  more  plane 
waves  in  our  group,  one  can  get  negative  refraction  of  the 
group  in  addition  to  the  energy.  Actually,  just  one  additional 
plane  wave  can  achieve  that.  Thus,  let  us  include  three  plane 
waves,  whose  wave  vectors  form  a  triangle,  rather  than  being 
parallel.  Let  the  magnitudes  of  the  wave  vectors  be  k,  k 
+  Skx ,  k+  $k2’  and  their  angles  with  the  normal  to  the  in¬ 
terface,  be  6 ,  6+  Sdi ,  0+  S02.  Inside  the  PIM  or  the  NIM, 
we  have 
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E  ye‘(k*x+kzz  "^{l +exp[/(«  ct)Sk\  +  ivkS6i] 

+  exp[((w  ct)8k2  +  ivk802]},  (14) 

with  u  x  sin  0+z  cos  0  and  v  x  cos  0  z  sin  0  for  the  PIM 
and  u  x  sin  0+az,  v  x  cos  0+bz,  and  kz  replaced  by  krz 
for  the  NIM.  Then  the  lines  whose  equations  are  u  constant 
and  v  constant  are  perpendicular  for  the  PIM.  Here  use  has 
been  made  of  the  following  expansion: 

krz(k  +  Sk,0+  80)^krz+a8k  +  bk80  (15) 

with 

a  k{m^0+c\nr\lvr)l\kr^, 
b  k  sin20/2|krz|. 

Here  vr  c(dnru>/dcti)  The  dependence  of  krz  on  k  and  0 
is  obtained  from  Eq.  (2).  The  condition  for  maximum  inten¬ 
sity  for  the  quantity  in  brackets,  the  long-wavelength  enve¬ 
lope  of  the  packet,  is  determined  by  the  equations 

(«  ct)Ski  +  vkS0i  2m\ir, 

( u  ct)8k2  +  vk802  2  m2TT, 

whose  solution  in  the  PIM  is 

x  (c i  sin0+c2  cos  0)  +  sin  0  ct, 

z  (c[  cos  0  c2  sin  0)  +  cos  0ct 

with 

Ci  2ir{m280\  8028k\), 

c2  2v(mi8k2  m28k\)/[k(80l8k2  8028k\ )], 

which  are  clearly  only  defined  for  8kx  !  8k2+  S0l  180 2. 

Inside  the  NIM,  the  solution  for  the  location  of  the  inten¬ 
sity  maxima  is 

x  ( c2a  C\b  bct)/(a  cos  0  b  sin  0), 

z  ( ci  cos  0  c2  sin  0+cos  0ct)/(a  cos  0  &sin0). 

From  the  expressions  for  a  and  b  under  Eq.  (15),  one  has 
a,b> 0  and  a  cos  0  b  sin  0>O.  Then  from  the  above  ex¬ 
pressions  of  x (t)  and  z(t),  one  has  dx/dt< 0  and  dzldt 
>0.  Thus  the  refraction  will  be  negative.  Let  the  angles  of 
the  line  u  constant  and  v  constant  in  the  NIM  with  the  z 
axis  be  a  and  >S,  respectively.  Then  one  has  tana 
a/sin  0,  and  tan /3  b! cos  0  kxtkn.  So  one  always  has 
7r/2<  a</3<  7r  inside  the  NIM.  From  the  above  expres¬ 
sions,  one  can  see  that  the  maxima  move  in  the  /?  direction, 
that  is,  antiparallel  to  kr .  The  velocity  of  interference  pattern 
in  the  NIM  defined  as  the  velocity  of  maximum  is  given  by 
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FIG.  3.  Electric  field  Re  E  of  negative  refraction  of  four  plane 
waves  with  wave- vector  magnitudes  k  8k,  k,  k+  8k,  k,  and  inci¬ 
dent  angles  0,  0  80,  0,  0+  80,  respectively.  Arrows  indicate  the 
directions  of  motion.  The  center  wave  number  is  k  \[5  with  inci¬ 
dent  angle  0  ir/6,  8k  0.2,  and  80  it/ 45.  Up  to  the  first-order 
approximation,  the  electric  field,  Poynting  vector,  and  the  momen¬ 
tum  of  this  group  of  plane  waves  are  Er  2e'<^(cos'nr-t-cos  8</>r), 
(Sr)  2 ( cos  nr + cos  8<j>^)\rlo),  P^e"  2Ak, /w,  respectively. 

This  velocity  is  independent  of  how  the  incident  wave  packet 
is  constructed.  The  refraction  of  a  group  constructed  from 
four  plane  waves  is  shown  in  Fig.  3  [8],  The  arguments  pre¬ 
sented  above  demonstrate  that  for  any  group  consisting  of 
three  or  more  plane  waves  whose  wave  vectors  are  not  col- 
linear,  the  group  refraction  is  negative. 

While  the  simulations  in  Fig.  3  clearly  show  that  the  in¬ 
tensity  maxima  refract  negatively,  the  normal  to  the  planes  in 
which  these  intensity  maxima  lie  are  directed  in  a  positive 
refraction  direction.  Thus,  if  one  was  to  imagine  smoothing 
out  all  intensity  variation  in  the  planes,  the  planes  would 
appear  to  refract  in  a  positive  direction.  We  believe  that  this 
is  a  remnant  of  the  positive  refraction  of  the  planes  of  inten¬ 
sity  maxima  [the  cosine  function  in  Eq.  (11)]  found  for  the 
interference  pattern  for  the  two  plane  waves  example  of  Ref. 
[3].  When  there  are  only  two  plane  waves,  this  is  the  only 
group  motion  that  we  see  in  the  NIM  since  for  a  group  con¬ 
sisting  of  two  plane  waves,  there  are  no  intensity  variations 
in  these  planes. 

Let  us  also  look  at  the  energy  flow  which  is  represented 
by  the  Poynting  vector.  For  three  wave  vectors  with  wave- 
vector  magnitudes  k  8k,  k  ,  k  8k,  and  the  angles  with  the 
normal  0  80,  0,  0+  80,  respectively,  the  magnetic  field 
can  also  be  calculated  from  Eq.  (14)  using  Maxwell’s  equa¬ 
tions  and  the  resulting  Poynting  vector  up  to  the  first  order  in 
both  8k  and  80  is  given  by 

(Sr)  ](l+4  cos2nr  +  4  costrr  cos  8rf>r)kr/(o 
k80  sinur  sin  <5^>,.(cos  0\+bz)/(o 
+  2 8k  cos2t u(a  krz/k) z/co,  (17) 


vr  vrkr . 


(16)  where 
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nr  (bz  + cos  Ox)k 80, 

8<j>r  (flz  +  sin  0x)8k  Scot. 

The  time  average  is  performed  over  the  period  corresponding 
to  the  average  frequency  of  the  three  plane  waves.  Here,  (Sr) 
is  not  localized;  rather  it  forms  a  lattice.  A  unit  cell  is  defined 
as  the  region  in  which  tit  changes  by  7T  and  8<f>r  changes  by 
27 t,  as  is  obvious  from  the  expression  for  Er  or  (Sr).  The 
area  for  each  unit  cell  in  NIM  is 

A  27T2/(acos0  b  sin  0)k8kS0. 

Instead  of  integrating  over  all  space  which  will  diverge,  one 
can  calculate  the  electromagnetic  momentum  for  each  cell. 
Ignoring  higher-order  terms  in  8k  and  80,  we  get 

Pf"  3A(  1  +2  8k/3k)Kr/2o)  (18) 

with  Kr  kr  2(x  sin  6+za)  8k/3,  the  average  of  the  three 
wave  vectors  which  make  up  the  group. 

A  packet  constructed  from  a  finite  number  of  plane  waves 
will  always  give  a  collection  of  propagating  wave  pulses 
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with  the  area  of  the  unit  cell  inversely  proportional  to  8k  and 
80.  For  the  above  localized  waves  made  of  finite  number  of 
plane  waves,  the  group  velocity  \gr  is  parallel  to  Pr  and 
antiparallel  to  the  average  wave  vector  Kr . 


IV.  CONCLUSION 

In  this  paper,  we  have  shown  that  for  any  localized  wave 
packet,  the  refraction  at  an  interface  between  a  PIM  and  a 
NIM  is  always  negative.  As  pointed  out  earlier,  it  is  essential 
for  a  correct  treatment  of  this  problem  to  use  wave  packets 
which  are  localized  in  all  directions  since  the  electromag¬ 
netic  field  from  any  physical  source  is  a  localized  wave 
packet. 
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Comment  on  “Wave  Refraction  in  Negative-Index  Materials:  Always  Positive  and 

Very  Inhomogeneous” 


Valanju,  Walser  and  Valanju  (VWV)  [1]  have  shown 
that  for  a  group  consisting  of  two  plane  waves  incident 
on  the  interface  between  a  material  of  positive  refractive 
index  (PIM)  and  material  of  negative  refractive  index 
(NIM),  the  group  velocity  refracts  positively.  Here  we 
show  that  this  is  true  only  for  the  special  two  plane  wave 
case  constructed  by  VWV,  but  for  generic  localized  wave 
packets,  the  group  refraction  is  generically  negative. 

The  sum  of  two  plane  waves  of  wavevector  and  fre¬ 
quency  (ki,uq)  and  (k2, W2),  considered  by  VWV,  can 
be  written  as  2e^k°  r_a'0^  cos[(l/2)(Ak  ■  r  —  Awt)], 
where  (ko,  wo)  the  average  wavevector  and  frequency  and 
(Ak,  Acj)  denote  their  differences.  Clearly,  the  argument 
of  the  cosine  is  constant  along  planes,  which  propagate 
in  time  along  the  direction  of  their  normal,  Ak.  We 
have  carried  out  numerical  simulations  of  wave  packets 
incident  on  the  PIM-NIM  interface  and  for  the  2-wave 
case  arrive  at  conclusions  similar  to  VWV.  For  arbitrary 
number  of  incident  plane  waves  whose  k  vectors  are  all 
parallel,  the  group  refraction  is  again  positive.  Note  that 
in  all  these  special  cases  the  packet  remains  nonzero  on 
infinite  planes. 

Here  we  show  that  for  any  wave  packet  that  is  spatially 
localized,  the  group  refraction  is  generically  negative. 

For  3  (or  more)  waves  whose  wave  vectors  not  aligned, 
the  group  refraction  will  be  negative.  For  example,  con¬ 
sider  three  wave  vectors  in  PIM  in  the  x- 2-plane,  whose 
magnitudes  are,  k  —  5k,  k,  k  —  5k  and  whose  angles  with 
the  2-axis  are,  6  —  59,  6,  6  +  56,  respectively.  The  disper¬ 
sion  k2  =  (u>2  —  uj2)(u>2  —  u>2)/c2(u>2  —  oj'q)  were  used  for 
NIM.  The  results  are  shown  in  Fig.  1.  The  wave  packet 
refracts  negatively,  in  obvious  contrast  to  VWV.  As  we 
have  seen,  two  plane  waves  result  in  a  wave  packet-like 
structure  which  is  constant  along  planes;  the  addition  of 
a  third  wave  breaks  the  planes  into  localized  wave  packets 
which  refract  negatively. 

A  packet  constructed  from  a  finite  number  of  plane 
waves  will  always  give  a  collection  of  propagating  wave 
pulses,  as  seen  in  Fig.  1.  A  wave  packet  localized  in  one 
region  of  space,  as  occurs  in  all  experimental  situations, 
can  only  be  constructed  from  a  continuous  distribution  of 
wavevectors.  Consider  such  a  wave  packet  incident  from 
outside  the  NIM,  E  =  E0  f  d2kf( k  -  k0)ei(k  r-w(k)t\ 
where  w(k)  =  ck.  If  /(k  —  ko)  drops  off  rapidly  as  k 


moves  away  from  ko,  w(k)  can  be  expanded  in  a  Taylor 
series  to  first  order  in  k  —  ko  to  a  good  approximation. 
This  gives,  E  =  Eoel(ko'r~u(ka')t g(r  —  ctko/ko),  where 
g{ R)  =  fd?kf( k  — k0)ei(k-k°)  R.  Inside  the  NIM,  k  and 
ko  in  the  argument  of  the  exponent  get  replaced  by  kr 
and  kro  which  are  related  to  k  and  ko  by  Snell’s  law. 
Then  the  wave  packet  once  it  enters  the  NIM  is  given  by 

Er  =  El0ei^krO'T~u,^krO^tgr(r  —  vgrt),  (1) 

where  <?r(R)  =  /  d2kf{ k  —  ko)e®R'Kk-k°)'Vkkrl.  Here 
kro  denotes  kr  evaluated  at  k  =  ko  and  vgr  =  Vkrw(kr) 
evaluated  at  kr  =  kro.  Thus,  the  refracted  wave  moves 
with  the  group  velocity  vgr.  Evaluation  of  Eq.  (1)  for 
a  Gaussian  wave  packet  shows  that  the  incident  packet 
gets  distorted  but  the  maximum  of  the  packet  moves  at 
vgr.  For  the  case  of  an  isotropic  medium,  considered  by 
VWV  [1],  the  group  velocity  is  anti-parallel  to  the  wave 
vector  in  the  medium.  Hence,  the  group  velocity  will  be 
refracted  the  same  way  as  the  wavevector  is,  contrary  to 
the  claim  of  VWV  [1], 

Thus  VWV’s  statement  that  the  “Group  Refraction  is 
always  positive”  is  true  only  for  the  very  special  (and 
rare)  wave  packets  constructed  by  them  and  is  incorrect 
for  more  general  wave  packets  that  are  spatially  localized. 
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FIG.  1:  Negative  refraction  of  3  plane  waves  with  k  —  6.32, 
5k  =  0.32,  0  =  7r/4,  50  =  n/60,  wo  =  2,  u)b  =  8,  ivp  —  10, 
and  c  =  1.  (Inset  above)  Wavevectors  for  the  3  plane  waves 
in  the  PIM  (left)  and  NIM  (right).  The  thick  arrows  indicate 
the  wave  packet  propagation  direction. 
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